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CHAPTER  I 


INTRODUCTION 

Background 

During  the  past  fifteen  years,  much  outside 
pressure — Congressional  and  taxpayer — has  been  placed  on 
the  Department  of  Defense  to  improve  the  management  of  its 
materiel  acquisition  system  as  well  as  other  resource 
expenditure  programs.  This  pressure,  combined  with  the 
aggressiveness  of  individuals  involved  in  the  materiel 
systems,  has  generated  many  innovative  techniques  to 
improve  the  support  of  our  weapons  systems  (1:22-23).  As 
a result,  the  focus  has  been  on  obtaining  the  maximum 
national  defense  for  minimum  dollar  expenditure;  a goal  to 
be  achieved  through  the  effective  use  of  national  resources. 

The  management  of  the  Air  Force  portion  of  these 
national  responses  is  the  responsibility  of  the  Air  Force 
Logistics  Command  (AFLC) , which  "...  must  maintain  the 
ability  to  replenish  base  stock  levels  through  a constant 
flow  of  recoverable  materiel  to  and  from  our  global  deployed 
Air  Force  units  [12:2-2]."  The  materiel  support  system 
required  to  accomplish  this  mission  involves  the  management 
of  1.8  million  line  items  worth  approximately  9 billion 
dollars  (7:34) . 
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Statement  of  the  Problem 


Forecasting  for  and  acquiring  spare  parts  is  an 
important  facet  of  the  materiel  support  system.  Forecasts 
are  used  to  develop  dependable  systems  and  to  formulate 
decisions  concerning  the  technical  and  economic  feasibility 
of  retaining  a system  in  use  (4:15-16).  The  economics  of 
spiraling  prices  and  austere  budgets  make  reliable  esti- 
mates, or  forecasts,  especially  critical.  The  exponential 
acceleration  of  spares  cost  provides  sufficient  justifica- 
tion for  basing  requirements  computations  upon  the  best 
available  forecasting  methodology  (1:22).  As  new  data 
sources  become  available,  as  improvements  are  made  in  fore- 
casting techniques,  as  added  insight  of  needs  is  received, 
and  as  other  uses  for  planning  factors  are  considered, 
current  forecasting  methods  require  modification  or  com- 
plete change.  Forecasting  is  an  evolutionary,  dynamic 
process,  and  the  method  which  served  well  yesterday  may 
not  continue  to  be  the  best  method  for  use  today  (8:351- 
358).  Operating  systems  managers  in  Headquarters  AFLC 
suspect  that  the  forecasting  technique  incorporated  in  the 
AFLC  materiel  support  system  is  biased,  i.e.,  the  expected 
value  of  the  forecast  error  is  not  zero.  The  problems 
analyzed  in  this  research  are  to  determine  the  bias  of 
the  current  AFLC  forecasting  technique  and  to  compare  this 
bias  to  that  of  other  selected  forecasting  methods. 
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AFLC  Forecasting  System 


The  Air  Force  Logistics  Command  manages  spares 
I according  to  the  provisions  of  the  Recoverable  Consumption 

Item  Requirements  System  (D041) , a computer  based  inventory 
control  system  (12:1-1).  In  this  context,  spares  are 
defined  as  recoverable  items  coded  with  an  Expendability- 
Recoverability-Reparability-Code  ( ERRC ) of  XD,  which  indi- 
cates that  the  spare  is  reparable  at  the  base  or  depot 
level  when  necessary.  Spares  constantly  move  through  a 
cycle  of  transactions  from  the  supply  activity's  warehouse, 
to  installation  on  an  end  item,  to  removal  from  a higher 
assembly  after  failure,  through  a repair  facility,  and  then 
back  to  the  supply  activity  ready  for  reissue  once  again. 
Spares  stockage,  after  initial  procurement  with  the  weapon 
system,  is  based  primarily  on  past  demands  or  reparable 
generations.  Information  concerning  these  reparable  gene- 
rations at  both  base  and  depot  level  is  accumulated  and 
consolidated  through  the  Air  Force  Recoverable  Asset  Manage- 
ment System.  This  data  is  input  to  the  D041  System  to 
provide  the  data  base  which  is  used  to  produce  the  fore- 
casts. These  forecast  computations  are  made  monthly,  for 
quarterly  periods.  The  initial  quarter  identifies  specific 
items  requiring  buy,  repair,  termination  and  disposal 
actions  (12:1-1).  The  second  quarterly  computation  (FY-2) 
is  used  for  adjusting  previously  initiated  item  actions, 
updating  workload  projections  for  annual  negotiation,  buy 
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actions  in  process,  termination  and  disposal  (12:1-1). 

The  third  quarterly  computation  (FY-3)  is  used  to  develop 
budget  requests  and  update  logistic  actions  in  process 
(12:1-2).  The  fourth  quarterly  computation  (FY-4)  is  used 
to  update  logistics  actions  in  process,  initiate  new 
actions,  and  provide  data  input  into  other  systems  (12: 

1-2)  . 

Forecasts  are  made  in  the  form  of  a reparable 
generation  rate,  which  provides  a basis  for  the  require- 
ments computation  performed  by  the  item  managers.  The  com- 
putations combine  the  forecast  with  other  asset  availabil- 
ity data  to  determine  actual  requirements.  These  require- 
ments are  further  modified  by  subjective  external  factors 
called  additives  (5). 

The  intricacy  and  subjectivity  of  the  requirements 
computation  procedures  are  beyond  the  scope  of  the  proposed 
research  effort.  This  research  focused  only  upon  the  sta- 
tistical model  incorporated  in  the  D041  System  used  to  fore- 
cast reparable  generations.  Consequently,  the  only  vari- 
ables considered  were  the  data  points,  or  actual  reparable 
generations  and  the  resulting  forecasts. 

Research  Justification 

- A number  of  sophisticated  statistical  forecasting 
techniques  are  available,  but  time  series  analysis  seems 
to  be  the  most  promising  for  forecasting  demand  data.  Time 


series  analysis  is  based  upon  the  concept  that  a sequence 
of  observations  composing  a time  series  is  one  realiza- 
tion of  jointly  distributed  random  variables.  The  past 
history  of  the  time  series  is  thus  used  to  provide  informa- 
tion about  the  mechanism  which  describes  its  evolution 
through  time,  and  the  mechanism  derived  is  used  to  fore- 
cast the  future  (1:18-19). 

The  procedural  forecasting  instructions  outlined 
in  current  Air  Force  directives  are  descriptive  transla- 
tions of  the  basic  mathematical  forecasting  models  used  to 
predict  future  requirements.  Little  management  attention 
is  given  to  the  forecasting  technique  unless  significant 
disparities  exist  between  forecasts  and  reparable  gene- 
rations. The  tendency,  even  then,  is  to  modify  the  fore- 
casts rather  than  to  question  the  validity  of  the  tech- 
nique generating  the  forecasts  (6) . The  effectiveness  of 
the  AFLC  spares  support  system  depends  upon  the  accuracy  of 
forecasts;  therefore,  it  is  imperative  that  the  most  accu- 
rate forecasting  techniques  available  be  used.  The  system 
currently  in  use  (D041)  depends  primarily  on  one  of  the 
simplest  statistical  techniques  available — the  weighted 
moving  average,  specifically  a twenty-four  month  moving 
average  (12).  Its  effectiveness  must  be  determined  and 
compared  to  that  of  other  forecasting  techniques  if  the 
Air  Force  is  to  operate  an  efficient  logistics  system. 
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The  moving  average  technique  smooths  fluctuations 


in  time  series  and  is  applied  to  trend  lines  as  well  as 
seasonal  and  cyclical  variations.  Smoothing  takes  place 
by  fitting  a straight  line  to  successive  groups  of  n years 
by  the  method  of  least  squares  (13:337).  However,  when 
n is  24  months,  a straight  line  may  not  adequately  fit. 

A parabola  may  perhaps  be  the  best  fit  in  the  least- 
squares  sense  (13:861).  The  moving  average  technique  also 
assumes  that  successive  observations  are  statistically 
independent,  i.e.,  assuming  this  year's  reparable  genera- 
tions are  independent  of  last  year's  generations,  but 
usually  this  is  not  true  (13:865).  This  technique  assumes 
that  the  data  are  stationary  (data  mean  and  variance  do 
not  change  over  time)  and  that  the  data  are  not  autocor- 
related  (independent  of  previous  time  periods) . The  moving 
average  technique  is  simple,  and  easy  to  use  because  of 
its  straightforward  calculations.  In  formula  form,  the 
predicted  value  of  reparable  generations,  g,  for  period  t 
(made  in  period  t-1)  is  such  that  the  average  for 

time  t+1  is: 

*t+l  = gt-l  + n (gt'^t-n)  (Equation  1) 

wherfe: 

A 

g^.  is  the  forecast  reparable  generations  for 
period  t+1, 


n is  the  number  of  periods, 
t is  the  rime  period. 


g is  the  actual  reparable  generation  in  period 
t,  and 

g.  is  the  actual  reparable  generations  in  period 
t-n. 

The  error  associated  with  this  forecast  is  the  difference 
between  the  forecast  value  and  the  actual  demand  value. 
With  the  simple  moving  average  forecasting  technique  the 
forecast  error  is  small,  but  its  standard  deviation  is 
greater  than  the  standard  deviation  of  the  actual  data 
(11:4).  The  moving  average  technique  merely  smoothes  out 
the  fluctuations  in  the  data  by  "moving"  the  arithmetic 
mean  value  through  the  time  series  (10:234).  This  tech- 
nique reacts  slowly  to  sudden  changes  in'  the  data  and  is 
ineffective  if  the  data  is  erratic.  If  a step  change 
in  demand  occurs,  the  moving  average  takes  n periods 
(the  moving  average  base)  to  level  off  at  the  new  demand 
(11:5).  However,  in  the  case  of  continuously  increasing 
demand,  estimated  demand  responds  slowly  and  always  under- 
estimates actual  demand.  When  this  occurs,  the  forecast 
becomes  more  biased  (11:4);  that  is,  the  absolute  value  of 
the  expected  value  of  the  error  increases.  The  forecast 
lags 'behind  the  trend  in  demand,  and  the  amount  of  lag 
depends  on  the  number  of  periods  in  the  forecast.  Altera- 
tions are  then  required  to  offset  the  bias  in  the  moving 
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average  forecast  causing  yet  another  problem.  The  "cor- 
rected" level  often  overshoots  the  actual  demand  and  does 
not  return  to  the  mean.  The  moving  average  technique  is 
not  easily  adapted  to  adjust  for  forecasting  over  long 
periods,  nor  does  it  allow  for  averaging  over  arbitrarily 
varying  time  periods  (11:6). 

Application  of  the  moving  average  technique  to 
time  series  data  requires  data  with  a characteristically 
linear  function  and  a definite  rhythmic  pattern  of  fluctua- 
tion which  repeats  periodically  (10:234). 

In  summary,  the  moving  average  is  simple  to  compute 
and  has  merit  if  the  basic  assumptions  have  been  met, 
but  it  has  a serious  disadvantage  in  that  it  responds  slowly 
to  changes  in  the  demand  pattern  (11:6).  Data  at  the 
beginning  of  a series  are  lost.  Moving  averages  may  gene- 
rate cycles  or  other  movements  which  are  not  present  in  the 
original  data.  Moving  averages  are  affected  by  extreme 
values  (11:6).  A method  of  averaging  that  allows  for  an 
easy  change  of  the  series  length  and  overcomes  the  effects 
of  extreme  values  should  be  used.  Further  support  for  dis- 
continuing the  use  of  the  simple  moving  average  was  found 
in  the  Predictive  Techniques  Study,  Phase  I:  "Comparison  of 
Some  Forecast  Methods"  done  by  M.  D.  Lum,  L.  L.  Blair  and 
J.  R.  Stuart.  They  used  generated  data  and  concluded  that 
the  24  month  single  item  moving  average  for  items  with 
long-term  exponential  or  linear  trend  was  the  worst  of  all 
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methods  tested,  and  that  double  exponential  smoothing  with 
specified  constants  using  Brown ' s tracking  signal  (an  out- 
of-control  signal  that  indicated  the  difference  between 
the  forecast  and  new  datum  has  become  significantly  larger 
than  one  would  expect  by  chance)  ranked  best  for  the  items 
with  exponential  trends.  The  study  also  recommended  that 
efforts  be  made  to  investigate  applications  of  other  tech- 
niques, and  that  real  data  be  used  in  the  future  studies. 

A more  reliable  method  for  forecasting  appears  to  be  neces- 
sary and  a search  for  this  method  has  led  statisticians 
to  the  exponential  smoothing  technique  (9:iii-v). 

Exponential  smoothing  uses  the  fundamental  idea 
that  any  new  estimate  of  demand  is  related  to  the  previous 
estimate  corrected  for  new  information.  Like  the  simple 
moving  average,  exponential  smoothing  also  assumes  station- 
aritv  and  independence.  The  estimating  equation  is: 

New  Estimate  = ©(New  Demand)  + (1-0)  (Old  Estimate) 

(Equation  2) 

where  0 < 0 < 1 is  the  smoothing  constant.  Small  values  of 
the  smoothing  constant  (0)  can  be  used  to  place  heavier 
emphasis  on  demands  from  the  distant  past  and  little  weight 
on  more  recent  demands.  Conversely,  larger  values  of 
theta  can  be  used  so  that  the  distant  past  has  little  effect 
on  the  forecast  (11:7).  By  varying  theta,  the  item's 
demand  history  can  be  increased  or  decreased  in  importance 
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for  making  forecasts,  and  the  method  can  be  repeated  using 
as  many  different  thetas  as  desired.  Overall,  when  the 


4 

smoothing  constant  is  small,  the  function  behaves  like  the 
average  of  a large  amount  of  past  data,  and  therefore  the 
variance  of  the  estimate  is  small.  When  the  smoothing 
constant  is  large  the  function  will  respond  rapidly  to 
changes  in  the  pattern.  The  advantages  of  exponential 
smoothing  are  simplicity,  ease  of  calculations,  minimal 
data  requirements  once  the  forecasting  method  has  been 
developed,  simple  control  procedures  for  the  forecasts; 
and  the  straightforward  manner  in  which  data  enter  the 
forecast  (11:8).  Therefore,  if  forecasting  is  performed 
several  times,  over  many  items,  exponential  smoothing  can 
be  routinely  accomplished.  For  unexpected  sudden  changes 
in  the  demand  pattern,  the  estimator  is  unbiased  in  the 
long  run;  but  for  steadily  increasing  demands  the  esti- 
mator is  biased  because  the  forecast  lags  behind  the 
demand  (11:9) . 

Double  exponential  smoothing  is  similar  to  single 

exponential  smoothing  and  the  conclusions  can  be  extended 

to  any  order  of  smoothing.  The  single  smoothed  value  of 

the  data  will  increase  to  the  value  of  theta  in  response 

to  the  data  input,  while  the  response  of  the  double 

2 

smoothed  value  is  theta  squared  (0  ) . The  order  of  smooth- 
ing is  based  on  the  choice  of  the  degree  of  the  polynomial 
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model  that  is  to  represent  the  data  and  the  subsequent 
selection  of  a smoothing  constant  for  the  respective  order. 

This  smoothing  constant  controls  the  number  of 
past  observations  that  have  any  effect  on  the  forecast. 
Flexibility  of  changing  the  constant  should  be  considered 
to  obtain  stable  forecasts  from  stable  processes  and 
forecasts  that  respond  rapidly  when  the  value  of  any 
coefficient  appears  to  have  changed  greatly  (11:9).  Two 
values  should  be  selected  for  the  smoothing  constant:  one 
for  stable  conditions  and  a higher  value  for  conditions 
of  rapid  change.  A small  value  includes  many  observations; 
but  for  larger  values,  only  recent  data  are  included. 

With  the  higher  order  models,  more  coefficients  must  be 
estimated.  Since  more  coefficients  are  needed,  the  number 
of  past  observations  should  also  be  increased  to  insure 
better  coefficient  estimators  by  using  as  many  data  points 
as  necessary  to  satisfy  the  degrees  of  freedom  requirement 
for  the  model. 

The  moving  average  method  does  not  attempt  to 
assess  rate-of-change  in  the  data  points  as  a function  of 
time.  In  other  words,  the  technique  only  estimates  the 
current  level  of  a factor  without  regard  for  whether  that 
factor  is  increasing  or  decreasing,  making  no  attempts  to 
adjust  for  trends.  The  forecasts  with  the  moving  average 
are  always  the  same,  regardless  of  the  forecast  lead  time 
(number  of  time  periods  into  the  future  for  which  the 
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current  forecast  is  being  made) . Thus  the  moving  average 
forecasts  straight  ahead.  In  contrast,  the  exponential 
smoothing  method  estimates  both  the  current  level  of  a factor 
and  the  current  linear  rate-of-change  for  the  same  order. 
Forecasts  made  are  a linear  function  of  the  forecast  lead 
time  using  these  estimates.  Exponential  smoothing  attempts 
to  account  for  and  to  follow  trends,  and  it  can  follow  the 
direction  of  these  trends  by  forecasting  upwards,  downwards 
or  straight  ahead  (9:9-10). 

Headquarters  AFLC  suggests  that  the  double  exponen- 
tial smoothing  method  should  be  analyzed  for  possible  AFLC 
use  (5)  and  that  the  present  D041  prediction  method  should 
also  be  compared  against  other  more  currently  developed 
predictive  techniques  (6) . 

This  thesis  continues  the  research  conducted  by 
Captain  Bruce  R.  Christensen  and  Mr.  Gene  J.  Schroeder,  who 
used  simulated  models  (2:Ch.l)  to  conclude  that  the  expected 
value  of  the  distribution  of  forecast  errors  presently  found 
in  the  D041  System  was  not  zero  (3:Ch.l-3).  Their  recom- 
mendation was  to  incorporate  actual  D041  data  into  the 
investigation  of  the  bias  produced  by  the  single  moving 
average  forecast  and  other  techniques  such  as  time  series 
analysis  (3:37) . 

Areas  of  Investigation 

There  are  two  areas  of  investigation  that  provide 
the  objectives  for  the  proposed  research: 
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1.  To  determine  if  the  D041  System  forecasts  are 

biased . 

2.  To  calculate,  review  and  compare  the  forecast 
error  of  time  series  analysis  methods  to  that  of  the  D041 
System  to  determine  which  of  the  techniques  is  the  best, 
in  terms  of  minimum  error,  to  forecast  failures  of  spare 
parts. 


Research  Question 

What  forecasting  technique  should  the  DC41  inventory 
management  system  use  to  provide  the  most  accurate  forecasts? 


Research  Hypotheses 

1.  The  use  of  the  weighted  moving  average  tech- 
nique leads  to  biased  forecasts  in  the  AFLC  material  sup- 
port system. 

2.  The  use  of  the  time  series  analysis  technique 
leads  to  forecasts  with  significantly  less  bias  than  that 
of  the  D041  System. 
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CHAPTER  II 


RESEARCH  METHODOLOGY 
Population  of  Interest 

This  research  focused  upon  the  characteristics  of 
the  D041  forecasting  technique,  the  technique  currently 
used  to  create  quarterly  forecasts  for  125,000  recoverable 
items  in  the  Air  Force  inventory.  Of  these  125,000  items, 
approximately  45,000  items  are  forecast  subjectively  by 
the  responsible  item  manager,  rather  than  by  the  D041 
System  based  on  past  demands  (5).  These  items  are  excluded 
from  consideration  since  they  do  not  typify  the  statistical 
technique  used  in  the  D041.  The  population  was  further 
limited  to  the  nine  Federal  Stock  Classes  listed  below: 


1560 

Aircraft  Structural 

1615 

Helicopter  Components 

1630 

Wheel  and  Brake  Components 

2915 

Engine  Components 

5841 

Airborne  Radar 

5865 

Electronic  Countermeasure 

5905 

Resistors 

5910 

Capacitors 

6615 

Airborne  Gyros 
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The  selection  of  these  classes  was  based  upon  the  past 
history  of  inaccurate  forecasts  and  consequential  problems 
as  recommended  by  Ms.  Anita  Clevenger  of  AFLC  (5).  Class 
1560,  Aircraft  Structural,  was  specifically  selected 
because  of  its  stable  demand  pattern.  The  population  was 
further  limited  by  the  requirement  for  eight  years  data 
necessary  to  examine  the  performance  of  the  techniques  over 
an  adequate  time  period. 

Definition  of  Variables 

Within  the  D041  forecasting  environment,  approxi- 
mately 20  variables  interact  to  influence  the  requirements 
computation  for  spares.  Examination  of  most  of  these 
variables  is  beyond  the  scope  of  this  research.  With  the 
focus  on  the  performance  of  the  forecasting  technique 
itself,  only  the  following  variables  were  considered  (5): 

1.  Base  reparable  generations,  i.e.,  the  number 
of  failures. 

2.  Organizational  Field  Maintenance  (OFM ) demand 
rate,  i.e.,  the  predicted  ratio  of  future  reparable  gene- 
rations to  the  predicted  number  of  flying  hours. 

3.  Operational  Program,  i.e.,  the  number  of  flying 
hours  for  the  spare  in  question. 

4.  Projected  Operational  Program,  i.e., 
of  flying  hours  forecast  for  future  periods. 


the  number 


Data  Requirements 


Data  used  in  this  research  consisted  of  31  quarters 
of  forecast  data  and  35  quarters  of  actual  reparable  gene- 
rations for  a single  Master  Stock  Number  (MSN) . To 
obtain  this  data,  three  elements  of  data  were  collected 
for  each  of  31  MSNs  over  the  time  period  FY  69-2  to  FY  77-4. 
These  data  elements,  as  shown  in  Table  2.1,  consist  of 
the  OFM  demand  rate,  the  base  reparable  generations  and  the 
past  operational  program.  To  obtain  a forecast  of  units 
demanded,  the  OFM  demand  rate  is  multiplied  by  the  pro- 
jected operational  program.  However,  the  forecasts  pro- 
vided by  AFLC  were  obtained  by  multiplying  the  OFM  demand 
rate  by  the  past  operational  program  (e.g.,  actual  hours 
flown).  This  data,  then,  does  not  reflect  the  actual  D041 
forecasts  and,  hence,  any  conclusions  drawn  concerning  the 
biasedness  of  these  forecasts  is  of  questionable  signifi- 
cance. 

Data  Collection 

The  required  data  was  available  from  the  D041  Tape 
Library  (5) . Computer  programs  had  to  be  written  to 
extract  the  required  information  from  the  D041  tapes  and 
transfer  it  to  a data  file  in  the  CREATE  System.  Appendix 
A describes  this  data  retrieval  in  detail.  Master  Stock 
Numbers  were  initially  selected  according  to  the  following 
edit  criteria: 
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1.  Element  is  in  one  of  the  9 classes  previously 


defined. 

i 

t 2.  Element  is  coded  XD. 

3.  Element  is  forecast  by  the  D041  Sytem,  not  by 

the  IM. 

The  initial  computer  run  resulted  in  500  MSNs  (Classes  1615, 
1630,  2915,  5841,  and  6615).  Subsequent  runs  against  the 
above  edit  criteria  reduced  the  usable  data  to  31  MSNs,  all 
in  class  6615. 

Research  Approach 

The  stated  purposes  of  this  research  are  to  analyze 
the  performances  of  different  techniques  used  for  fore- 
casting reparable  generations  and  to  compare  their  accurate- 
ness. The  basic  approach  was:  (1)  to  compare  the  D041 

single  moving  average  forecast  with  actual  reparable  genera- 
tions for  the  same  time  period  to  determine  the  forecast 
bias,  and  (2)  to  compare  the  two  methods  (D041,  and  time 
series  analysis)  with  each  other  to  determine  which  tech- 
nique offered  the  least  amount  of  bias. 

Testing  the  First  Research  Hypothesis 
The  first  research  hypothesis  was:  The  use  of  the 
weighted  moving  average  technique  leads  to  biased  forecasts 

in  the  AFLC  materiel  support  system.  This  is  to  say  that 
the  expected  value  of  the  forecast  error  presently  found 

18 


- '■*  i • »■**  % . yV 


4 


9 «•  ]<~*- 


I 


in  the  D041  System  is  significantly  different  than  zero. 
The  hypothesis  may  then  be  stated  as: 

Hq  : E(eT)  = 0 

H,  : E(e_)  * 0 

± L 

where: 

E(eT)  is  the  expected  value  of  the  forecast  error, 
for  fixed  lead  time  r • 

Each  D041  forecast  from  quarter  five  (FY  70-2)  to  quarter 
34  (FY  77-3)  was  compared  to  the  actual  reparable  genera- 
tions for  quarter  six  (FY  70-3)  to  quarter  35  (FY  77-4). 
The  forecast  error  associated  with  each  time  period,  for 
each  Master  Stock  Number  was  calculated  as: 

eT(t)  = g^(t)  - g ( t ) (Equation  5) 

where : 

e (t)  is  the  error  term  associated  with  (t) , a 
lead  time  of  T, 

g (t)  is  the  forecast  for  period  t made  t periods 
earlier,  and 

g(t)  is  the  actual  reparable  generations  in 
period  t. 

The  final  data  pattern  consisted  of  four  error  distribu- 
tions (one  for  lead  times  of  1,  2,  3,  and  4 quarters)  for 
each  'Master  Stock  Number.  The  data  points  were  the  error 
[ e _ ( t ) ] values  for  t = 1,  2,  3,  ...  x for  a fixed  lead 


Table  2.2  should  clarify  the  explanation  of  fore- 
casts made  for  various  time  periods  for  lead  times  one  and 
two  (3:20) . 

The  first  research  hypothesis  was  designed  to  deter- 
mine if  the  D041  forecasting  method  is  biased.  The  expected 
value  of  each  error  distribution  for  a fixed  lead  time  was 
calculated  according  to 

n 

E(ex)  = [g_(t)  - g(t)]  (Equation  6) 

t=l 


where : 

E(eT)  is  the  expected  value  of  the  forecast  error 
made  for  lead  time  t, 

g (t)  is  the  forecast  reparable  generations  for 
period  t with  lead  time  t, 

g(t)  is  the  actual  reparable  generation  in  period 
t , and 

n is  the  number  of  data  points  in  the  distri- 
bution. 

The  hypothesis  was  then  tested  for  each  error  distribution. 
Rejection  of  the  null  hypothesis  means  that  the  D041  fore- 
cast technique  is  biased  (13) . 

A computer  program  was  written  to  compute  the  fore- 
cast errors,  and  an  SL  library  program  (SIMFIT)  was  used  to 
determine  the  statistic  properties  of  these  errors. 


20 


4 


TABLE 


Testing  the  Second  Research  Hypothesis 


The  second  hypothesis  states:  Use  of  the  time 
series  analysis  techniques  leads  to  forecasts  with  signifi- 
cantly less  bias  than  that  of  the  D041  System.  To  obtain 
a basis  for  comparison,  forecasts  were  made  from  the  same 
data  base  and  lead  times  using  SL  library  program  TCAST. 
These  forecasts  were  compared  to  the  actual  reparable 
generations  for  the  period  and  the  resulting  errors  were 
calculated  according  to  Equation  5.  The  expected  value  for 
each  distribution  was  then  calculated  according  to  Equation 
6 and  tested  with  the  following  hypothesis: 


H,  : E(e_)  ? 0 

1 i 

As  with  the  first  hypothesis,  SL  library  program 
SIMFIT  was  used  to  compute  the  statistical  properties  of 
these  errors.  If  the  statistical  test  resulted  in  the 
rejection  of  the  null  hypothesis,  it  was  inferred  that  the 
tested  forecasting  technique  was  biased.  Conversely,  if 
the  null  hypothesis  could  not  be  rejected,  the  inference 
was  that  the  tested  forecasting  technique  is  not  biased. 

The  expected  values  of  the  error  distributions 
obtained  from  the  two  forecasting  techniques  were  tested 
for  significant  statistical  difference  according  to  the 
following  hypothesis: 
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If  the  statistical  test  results  in  the  rejection  of  the 
null  hypothesis,  the  time  series  (TCAST)  forecasting  tech- 
nique is  inferred  to  be  more  accurate  than  the  D041  tech- 
nique. Conversely,  if  the  null  hypothesis  cannot  be 
rejected,  the  inference  is  that  there  is  no  significant 
difference  between  the  techniques. 

List  of  Assumptions 

1.  The  D041  forecasting  procedures  have  been 
applied  consistently  through  the  period  of  time  examined 
by  this  research. 

2.  The  error  terms  for  a specific  lead  time  are 
statistically  independent. 

3.  Forecasting  errors  are  normally  distributed. 

List  of  Limitations 

The  research  was  limited  to  the  forecasting  tech- 
niques as  applied  to  a subset  of  reparable  items.  Extrapo- 
lation to  forecasting  for  other  non-reparable  items  may 
require  additional  research. 
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CHAPTER  III 


FINDINGS  AND  ANALYSIS 

The  D041  Forecasting  Method 
The  D041  forecasting  method  was  investigated  to 
determine  if  the  forecast  was  biased.  The  forecast  error 
distributions  and  their  statistics  were  computed  by  com- 
puter programs  for  lead  times  of  one,  two,  three,  and  four 
quarters  (refer  to  Table  D.l).  For  this  particular  MSN 
(661500735893) , the  null  hypothesis  could  not  be  rejected 
at  any  reasonable  confidence  level  for  any  of  the  four  lead 
times  examined,  inferring  that  the  D041  forecasting  tech- 
nique is  unbiased  for  these  data. 

The  graph  plots  (Figures  C.9  through  C.12)  show  the 
error  distribution  to  approximate  rather  flat  normal  curves 
which  are  in  fact  very  close  to  uniform  distribution  curves. 

The  TCAST  Forecasting  Method 
The  exponential  smoothing  forecasting  method  (TCAST) 
was  investigated  to  determine  if  it  produced  forecasts  which 
were  biased.  The  distribution  of  forecast  errors  and  its 
statistics  were  computed  for  lead  times  of  one,  two,  three, 
and  four  (refer  to  Table  D.2). 

The  TCAST  model  was  optimized  for  each  forecast 
made.  Nearly  72  percent  of  the  102  forecasts  made  were 
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optimized  with  third  order  exponential  smoothing,  which 
takes  into  account  the  quadratic  rate  of  change  in  the 
single  exponentially  smoothed  average-  Third  order  exponen- 
tial smoothing  normally  results  in  more  responsive  fore- 
casting than  first  or  second  order. 

Similarly,  larger  alpha  (smoothing  constant)  values 
usually  result  in  more  responsive  forecasting,  but  in  these 
data,  86  percent  of  the  102  forecasts  were  optimized  with 
smoothing  constants  of  . 1 or  .2  (refer  to  Table  3.1). 

The  null  hypothesis  could  not  be  rejected  at  any 
reasonable  level  of  confidence  for  any  of  the  lead  times 
of  one,  two,  or  three,  inferring  that  the  TCAST  method  of 
forecasting  is  unbiased  for  these  data.  However,  for  a lead 
time  of  four  quarters,  the  null  hypothesis  could  be  rejected 
with  87  percent  confidence,  inferring  that  the  TCAST  method 
of  forecasting,  for  a lead  time  of  four  quarters,  is  biased 
for  these  data. 

The  graph  plots  (Figures  C.13  through  C.16)  of  the 
error  distributions  approximate  rather  flat  normal  curves 
which  are,  in  fact,  very  close  to  uniform  distribution 
curves . 


Comparison  of  the  Two  Forecasting  Methods 
' The  statistics  of  the  D041  and  the  TCAST  forecast- 
ing methods  were  compared  to  determine  if  the  TCAST  method 
offered  less  bias  than  the  D041.  The  null  hypothesis  could 
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not  be  rejected  with  any  reasonable  confidence  for  any  of 
the  four  lead  times,  inferring  that  there  is  no  statistically 
significant  difference  in  the  expected  values  of  the  fore- 
cast errors  for  these  two  methods. 

For  all  four  lead  times,  the  D041  forecasting  method 
produced  positive  mean  error  values,  inferring  a somewhat 
consistent  tendency  to  over-forecast.  However,  the  TCAST 
method  produced  negative  forecast  error  means  for  lead  times 
one,  three  and  four,  inferring  a tendency  to  under-forecast. 


CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

The  first  research  hypothesis  stated  that  the  use 
of  the  weighted  moving  average  forecasting  technique  by 
the  AFLC  material  support  system  led  to  biased  forecasts. 
Because  of  the  nature  of  the  forecasts  provided  by  AFLC, 
no  conclusions  can  be  reached  concerning  the  biasedness 
of  the  D041  forecasting  technique. 

The  second  research  hypothesis  stated  that  the  use 
of  time  series  analysis  technique  would  lead  to  signifi- 
cantly more  accurate  forecasts.  Again,  no  firm  conclu- 
sions can  be  stated  regarding  differences  in  accuracy 
between  the  two  techniques.  However,  the  TCAST  forecasts 
were  statistically  unbiased  and  made  no  use  of  "after  the 
fact"  information  as  did  the  D041  forecasts.  Also,  the 
variance  of  the  distribution  of  forecast  errors  for  the 
TCAST  technique  was  substantially  less  than  the  variance 
associated  with  the  distribution  of  forecast  errors  for 
D041.  This  smaller  variance  means  that  the  risk  associated 
with  'using  the  TCAST  forecast  is  less  than  that  associated 
with  using  the  DO 41  method. 
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Recommendations 


Recommendations  for  further  research  include  the 
following : 

1.  Research  should  be  conducted  using  a larger 
variety  of  items.  The  data  needed  for  this  research  is 
available  in  AFIT/CREATE  Library  under  permanent  D041 
file. 

2.  The  future  program  forecasting  procedures  and 
their  ramifications  should  be  studied.  The  projected  fly- 
ing hours  (future  program)  are  a driving  force  in  the 
calculations  to  determine  how  many  reparable  assets  of  a 
given  type  will  be  required. 

3.  The  Box  and  Jenkins  forecasting  technique  should 
be  investigated  as  a possible  better  forecasting  method. 

This  method  should  be  compared  to  both  the  D041  results  and 
the  exponential  smoothing  time  series  results  in  this 
thesis . 

4.  The  thesis  results  obtained  from  generated 
data  by  Captain  Eruce  Christensen  and  Mr.  Gene  Schroeder 
should  be  compared  to  the  results  of  this  thesis  and  any 
other  techniques  studied  using  real  data. 


APPENDIX  A 

D041  DATA  BANK  DATA  F.ETRIEVAL 
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The  Current  D041  Data  Bank  consists  of  approxi- 
mately 330  magnetic  tapes  containing  a staggering  volume 
of  historical  information.  It  is  in  fact  this  volume 
which  greatly  complicates  actual  usage  of  the  available 
data.  In  order  to  ferret  out  the  bits  of  information  per- 
tinent to  the  issue  at  hand,  hours  of  computer  time  were 
consumed  in  a computer  search  through  the  information. 

In  acquiring  the  data  for  this  research  effort, 
much  time  was  dedicated  to  locating  a knowledgeable  person 
to  explain  system  access  procedures.  Ms.  Anita  Clevenger 
of  the  Policy  and  Systems  Division,  HQ  AFLC  and  Ms. 

Virginia  Thomas  of  the  CREATE  Application  Support  Group 
were  able  to  jointly  and  cooperatively  furnish  the  required 
guidance.  Ms.  Thomas  can  furnish  copies  of  the  quarterly 
D041  Data  Bank  control  ledger  and  the  file  format  for  each 
record . 

The  first  set  of  problems  encountered  involving 
learning  the  use  of  tapes.  A set  of  permanent  tapes  was 
secured  through  the  use  of  System  NEOTAPE.  In  order  to 
preclude  confusion  and  the  inadvertent  destruction  (write- 
over)  of  one  of  the  data  tapes,  a control  log  was  used. 

For  each  program  run  the  numbers  and  labels  of  the  tapes 
read  from,  and  written  to,  were  recorded  in  this  control 
log.  A CARDIN  NOSLEW  control  card  was  used  immediately 
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prior  to  defining  the  "write  to"  tape  in  order  to  suppress 
any  carriage  control  format  instruction  and  therefore  to 
obtain  accurate  "read"  and  "write"  actions.  Some  of  the 
D041  tapes  are  in  marginal  condition,  which  occasionally 
caused  a run  to  abort  with  an  irrecoverable  input  error. 

In  each  case,  running  the  same  program  again  after  request- 
ing special  assistance  of  the  operator  resulted  in  a satis- 
factory run.  Should  this  have  failed,  Ms.  Thomas  could 
have  requested  the  reconstruction  of  the  defective  tape. 

Program  design  incorporated  the  minimum  number  of 
tapes  for  each  single  run  because  of  the  limited  number  of 
tape  drives  (14)  available  on  the  CREATE  System.  If  a 
single  run  required  more  than  three  tape  drives,  some  delay 
was  experienced.  Several  smaller  runs  rather  than  one 
large  one  were  preferred. 

In  this  research  effort,  the  basic  programming 
approach  was  to  run  the  program  against  the  first  available 
data  to  arbitrarily  select  900  Master  Stock  Numbers  (MSNs) 
according  to  certain  criteria,  and  write  them  to  a data 
tape.  The  next  run  read  these  MSNs  from  the  data  tape, 
read  through  the  next  quarter's  data  tapes  until  the  MSNs 
and  desired  list  of  information  matched,  and  then  wrote 
all  of  that  information  to  another  data  tape.  The  follow- 
ing run  read  that  information,  read  the  next  quarter's 
information,  and  wrote  all  the  desired  information  to  a 
third  data  tape.  In  other  words,  all  the  desired 
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information  was  transferred  back  and  forth  through  each 
consecutive  run.  From  this,  it  is  obvious  that  runs 
had  to  be  made  in  sequence.  From  discussions  with  the 
computer  operator,  it  was  learned  that  a CARDIN  MSG1  con- 
trol card  could  be  used  to  instruct  the  operator  to  run 
that  program  only  after  the  completion  of  a previous  run. 
Therefore,  an  entire  series  of  runs  could  be  entered  at 
one  time.  Preparing  the  runs  in  sequence  led  to  the  fact 
that  prior  to  quarter  74-3,  only  first  and  third  quarter 
data  was  retained  in  the  Data  Bank.  After  74-3,  data  was 
retained  for  each  quarter. 

Beginning  with  quarter  75-4,  MSNs  became  National 
Stock  Numbers  (NSNs)  including  the  "00"  (e.g.,  6615003147504). 

Prior  to  75-4  MSNs  were  really  Federal  Stock  Numbers 
(FSNs ) without  the  00  (e.g.,  66153147504) . In  order  to 
make  this  transition  in  the  series  of  runs,  a conversion 
program  was  written  to  read  the  information  already  col- 
lected from  a data  bank  and  write  it  back  to  another  data 
tape,  inserting  the  00  with  a format  statement.  See 
Figure  A.l  for  a copy  of  this  program. 

In  order  to  retrieve  the  data  necessary  for  this 
research,  three  specific  pieces  of  information  were  required, 
for  each  quarter,  for  an  unknown  number  of  MSNs.  The  OFM 
demand  rate  is  available  on  the  01  record;  the  base 

reparable  generations  are  available  from  the  03  record; 
and  the  past  program  is  available  from  the  11  record. 


The  first  run  selected  the  initial  set  of  (in  this  case) 
Federal  Stock  Numbers:  100  in  each  of  five  classes.  The 
only  other  criteria  were  a requirement  for  existing  data  in 
the  first  quarter  (70-2)  and  a computer  made  forecast  for 
that  quarter  (card  column  090  of  the  01  record  should  be 
blank) . An  "E"  in  that  card  column  would  indicate  a 
manipulated  forecast.  The  program  (Figure  A. 2)  was  run 
against  the  01  record  in  quarter  70-2,  which  consisted  of 
three  tapes.  These  tapes  appeared  to  be  in  stock  number 
sequence.  Results  were  500  Federal  Stock  Numbers  and  the 
associated  OFM  demand  rate  for  that  quarter.  This  informa- 
tion was  then  written  to  a tape  which  in  turn  provided  the 
basis  for  the  run  against  the  next  quarter's  data. 

These  runs  against  quarter  71-1  and  subsequent 
quarters  were  made  with  the  program  in  Figure  A. 3,  which 
is  the  final  run  against  quarter  77-4.  This  series  of 
runs  retrieved  the  OFM  demand  rate  for  each  of  the  stock 
numbers  remaining  in  the  data  sample.  Since  stock  numbers 
were  eliminated  when  no  data  was  found  or  when  a manual 
forecast  was  made,  the  number  of  stock  numbers  in  the 
sample  rapidly  decreased  from  500  to  46  with  the  final  run. 

With  the  demand  rate  data  in  hand,  the  next  runs 
were  made  against  the  03  and  11  records  to  obtain  the  past 
program  and  reparable  generations  data.  The  records  show- 
in  the  D041  Data  Eank  Control  Ledger  as  only  the  03  record, 
but  the  03  through  11  records  are  carried  on  those  tapes. 
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Sequence  appears  to  be  by  Air  Logistics  Center,  then  stock 
number,  then  record  type.  Since  eight  quarters  of  data 
are  carried  on  each  quarter's  records,  it  was  necessary 
only  to  run  against  the  71-1,  the  72-3,  the  74-3,  the  76-3 
and  the  77-4  tapes  to  obtain  a full  set  of  data  for  the 
stock  numbers  in  the  data  sample.  A sample  of  the  program 
can  be  seen  in  Figure  A- 4. 

Programming  logic  for  this  series  was  similar  to 
that  used  for  obtaining  the  OFM  demand  data.  The  end 
results  were  22  quarters  of  OFM  demand  data  (because  of 
the  missing  quarters)  on  one  tape  and  35  quarters  of  past 
program  and  reparable  generation  data  on  the  second  tape. 
To  mesh  them  together,  compute  the  forecast  and  write  all 
the  data  to  a permanent  data  file,  another  program  was 
written  (Figure  A. 5).  As  the  data  was  read  from  the  tapes 
prior  to  adjusting,  the  quarters  did  not  match.  The  data 
for  a sample  MSN  appeared  as  in  Table  A.l.  The  missing 
OFM  demand  rates  were  filled  in  by  straight  line  averaging 
the  nearest  actual  demand  rates.  The  limitations  of  this 
step  were  recognized,  but  it  was  deemed  more  valuable  to 
compute  the  statistics  for  a larger  sample  of  error  points 
than  to  rely  on  strictly  actual  available  data.  Since  no 
demand  rates  were  available  for  the  first  four  quarters 
(FY  69-2  through  FY  70-1)  no  forecasts  were  analyzed  for 
those  quarters. 
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The  end  result  of  this  program  run  was  a usable 
stream  of  35  quarters  of  data.  The  full  data  file  for  all 
31  remaining  stock  numbers  is  established  as  a permanent 
file  under  the  name  D041.  This  file  is  accessible  for 
other  potential  uses.  The  data  for  the  stock  number 
examined  in  this  research  is  tabulated  in  Table  A. 2. 

The  actual  data  points  to  be  examined  in  this 
research  were  the  errors  of  the  forecast.  A computer  pro- 
gram (Figure  A. 6)  to  read  the  D041  data  (forecast  and 
reparable  generations)  and  compute  the  forecast  errors 
for  lead  times  one,  two,  three,  and  four  (Table  A. 3)  was 
accordingly  written. 

Time  became  one  of  the  largest  problems  in  obtain- 
ing this  data  file.  The  programs  themselves  normally 
required  15  to  35  minutes  of  computer  time  each.  When  this 
requirement  is  combined  with  the  tape  drive  requirements 
and  sequencing  requirements,  a data  retrieval  project 
becomes  a formidable  task  requiring  a significant  number 
of  calendar  days  to  complete.  The  situation  is  complicated 
by  computer  down-time  and  the  fact  that  jobs  this  large 
are  only  run  at  night  and  on  weekends.  Once  the  tape  and 
formats  were  identified,  writing  and  running  the  set  of 
programs  in  this  appendix  required  approximately  six  weeks 
of  continuous  daily  runs. 
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1800  PRINT  I 25  3,  FSCC I,J)  ,FIIN( I,  J) ,OF"DR(  I,J) 

1613  I 250  FORMAT  ( 1 X , 1 4, 1 X , A 1 i , 2X,F6. 4 ) 

1320  NRITE  ( I 6,  I 255,  END=I  230)  FSCC  I , J)  ,F  I IHC  I , J)  ,OF?‘DR  ( I , J ) 
1 8 30  I 255  FORMAT  ( 14  , A 1 1 , F6.  4) 

1843  1257  CONTINUE 
1853  1260  CONTINUE 
I860  I 2/0  CONTINUE 
1873  12  30  CONTINUE 
1 3300 
I £9-3  STOP 
1903  END 
I 91  OS  J EXECUTE 
I 920$: LIMITS: 4 0 

1 933$ : TAPE : 1 3, XI D, , 70460 ,, 73-2  DD301 
1 943$: TAPE t I 4 , X2D , , 73363 , , 70-2  DDB01 
I 950$  J FF ILE» I 6,N0SLEN 
1 960$ : TAPE : I 6, A I D , , 7367 1 , , MSNDATA/R ING 
1970$: END JO 3 
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I •00ri  C FIGURE  A. 3 OF”  DEMAND  DATA  RHTR  I :V  AL  PROGRA-" 

1 . I 0 C 
I D20C 

1 03"  YVD , r?  ( SL ) 

1 ^40S:  IJENTs  .:P1  I 9^,  AFITSL  77 A BPAiNTLEY  LORE  4 Ail 

I 05:1$  * ‘(SGI  t'tUST  3E  PUD  AFTER  SNJM3  G5' 16  3 PLEASE 
I 060  S : DPT  IO:  i : FORTR  AN  , . ! O' DAP 
I 070$  s FORTY : .'IFOR  !■' . . i L-  IO 
I 030  C 

1 090C  RROGRA,"  DEMAND 
I I 00  C 

1110 C THIS  PROGRAM  SEARCHES  THROUGH  EACH  QUARTER'S  DATA 
II 20C  TO  RECORD  THE  OF”  DEMAND  RATE. 

I I 30  0 

1140  CHARACTER  FI  IN*  I I (54),FN*I  I fDC*l 
1150  INTEGER  FSC( 54 ) , CLASS,  Q,GF, OTR 
I 1 6"  DIMENSION  OFMDR(52,  22) 

1 1 70C 

I 1 3 0 C READ  IN  THE  VALUES  FDR  THE  TOTAL  DUMBER  Or  MSil 
I 1 90 C IN  EACH  CLASS. 

1 2 ODD 

1210  READ  3 00, D:JM 
1220  300  FORMAT  (V) 
i 2 3 * C 

1 24 CC  SET  THE  QUARTER  DUMBER  BEFORE  EACH  RUN . 

I250C 
1250  0=22 
I 270  C 

1 230 C INITIALIZE  THE  TAPE  DEVICE  TO  3E  READ  FROM  ID 
I290C  THE  DD4 1 DATA  BANK. 


310  ID 

T=  1 3 

*3  9 O f ' 

3 3 DC 

SET  TDE  VALUE  FOR  THE 

VARIABLE  REPRESENTING  THE 

340  C 

LAST  TAPE  DEVICE  ID  Tl 

JE  DATA  BANK  TO  bE  READ  FROM 

35  OC 

ON  Til S RUN. 

360  C 

37^  IE 

D=  1 3 

3 8 o : 

°9r'C 

INITIALIZE  THE  DEMAND 

DATA  ADRAYH  TO  -1.  ID  ORDER 

4 fWC 

TO  PROVIDE  INDICATORS 

FOR  DETER.M I. < IDG  .HI  CH  MSN 

4 I ‘'1C 

WERE  NOT  MATCHED. 

4R0C 

1 4 30  DO  I 120  J=l , NUV 
I 44  ' OFMDR  ( J,D)=-1  . 

I 450  C 

1 4500  READ  TUE  MSN  AND  DEMAND  DATA  FROM  THE  DATA  TAPE. 
I ■*■  7 1 3 


4 3'! 

RE 

AD 

( 1 6 

, 1 v 

t 

END=  1 

°2  *)  F SC  ( J ) , r 1 1 D ( J ) 

49/ 

1 .3 

.0  0 

FOR 

NAT 

( i 

4, All 

) 

50 

DO 

1 

09  3 

QTR= 

i , 

C-l 

51  " 

RE 

AD 

( 1 6 

*101 

• t 

END=1 

22  3)  OF  i.jR(  J ,QTR  ) 

520 

1 J 

1 0 

F )R 

MAT 

(F 

) .4  ) 

530 

1 3 

9.3 

COD" 

r i d u 

^2 
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o o o -o  o o o on  ( :i  ca  :ij  u.  rr>  cr  oj 


I 5b-  C 

I 56  ? C >2 ARCH  T-»RU  rHI5  QUARTERNS  DATA  AID  READ  THE  FS C 
157  k;  fun,  oe"  h.T'aijd  fate,  and  demand  ate. 

I 5 6 ' C 

15V  DO  1210  I I JX=  I , i EDO  DO 
15'*.'  GO  TO  I 1 29 
I 6 >00 

I620C  DETER*' I ME  IE  THE  LAST  AVAIL  AULA  TARE  HAS  SEEM  EE 


64.'  1125  I.i:=II!T+l 
>50  IE  ( INT.GT.  IH.I)  GO  fO  I "2  ' 

•56.*  1129  READ  ( IN T , 1 I 30 , EH J= I l 25 ) CLASS, r- , JR  , DC 
6 70  1 i 30  FI)  2 V AT  ( -> X , 1 4 , A I I , 65  X , F 5 . .. , A I ) 

630C 

690 C MATCH  T'E  ITEM  READ  TO  THE  oASIC  DATA  LIST 
700C  PREVIOUSLY  SAVED.  IE  THE  STOCK  NJM3HRS  HATCH, 
7IUC  I.1CRE."H.MT  THE  COTTER  (KN7R ) TO  TRACK  THE  .JLM3HR 
72 OC  OF  DATA  POINTS  AH  J RECORD  TUE  L"\iD  DATA.  It-  T 
7 30 C STOCK  MUM 32 i READ  DOESN'T  AT CH  ANY  I..  THE  DATA 
740C  LIST,  READ  ANOTHER  MJ MS HR . 

7 5 0 C 

760  IF  (CLASS.  IE.  6615)  JO  TO  1210 
770  DO  1135  J = !,Nj'< 

780  IF  ( F.< . ME  .F  1 1 ( J ) ) JO  TO  1165 
790  KMTR=KNT R+l 


HOC  IF  THE  DEMAND  CODE  IS  E,  DO  MOT  SA 
20  C OTHER'S  1 32,  SAVE  THE  DH”  AND  DATA. 

30  C 

40  IF  ( DC. ME.  ' 1*  ) JO  TO  I i 63 
5.-'  DO  TO  1214 

63  No'1  OF  MDR  ( J , D ) =D  R/  1 OR  "50  . 

70  IF  (KN’TR  .LT.NlI")  GO  TO  121" 

-30  GO  TO  1 22  0 
PC  I 1 65  CONTINUE 

10  1210  CONTINUE 


30C  UPDATE  THE  FILE.  FIND  EAC--  "SM  .IN"!  EE  'AND  DATA 
40 C STILL  EQUAL  TO  -i  . AMD  ELIMINATE  IT  ;r  ADJUSTING 
ore  THE  FIINS  FRO  V THAT  POI  if  ON.  T !E  : CHORE  N:  NT  7" 
60C  iJU?*BHR  Or  MS  M'S  R2MAI  ■ IN 3 IN  THAT  CLAES  AMD  AQJU 
70 C ALL  THE  DEMAND  DATA  COLLECTED  THUS  EAR  1 0 MAIN  7 A 


930  C 

ALIGN "ENT  . 

99  -9  C 

<r<  *• 

PRINT  THE  DAT 

■IOC 

C 20  1 

220  IT?=NJ” 

030  D 

•5  1 250  J = 1 , IIP 

204.:  PRINT  I 23  0,  ESC  ( J)  , E 1 1 M(  J ) 

■>:  5"  1 23J  FORM.  AT  ( 1 X,  14,  I a,  A 1 1 ) 

..RITE  (17.1  24  3.  END*  127*)  ES  C(  J ) . E 1 1. . ( J ) 


II 


2 'Q0  JO  I 24  3 G'i'-'=  1 , ! 

2G9  1 P3I.iT  I 245,OF!JRCJ,  OTP) 

2 1 20  1245  FOP  VAT  ( 1 1 X , F i> . 4 ) 

2112  • WITH  (17,1  247, HA  )=l  27  J)  OF  VJR ( J , JTP ) 

2120  124  7 FOR  VAT  ( Fo  .4  ) 

2130  124  3 COimWUE 
2140  GO  TO  1250 
2150  1250  COiiTI.NUH 
2160  PRINT  1260,  NGN 

2 1 7 0 C 1260  FORMAT  ( //  1 X , 4 I HTH  £ u'U  Vs  HR  OF  STOCK  liGMEHRS 
2I3'H  RH MAIN  IMG  IS  ,12) 

2190  1270  COMTI.NUH 

**  o , as* 

J* 

2210  STO D 

2220  h:o 

2230s  s HXEC'JTH 
2240$  sL  I M ITS  * 3 0,  I 5'C 
2250 3 : JATA ! 1+ 

2260  51 

22  703  sTAPH s I 3, X 1 D , , 735  I 2 , , 77-4  JD50 1 
22303 sTAPH s 1 6 , X20 , , 737  39  , , JATAF ILH7 
2290$  sFFILHs 1 7, MOSLEM 

2300$  sTAPHs  I 7,  A I C, , 7 J846  , , JATAF  I LH3/R  Ir.G 
23 1 0$  s HMDJOd 
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I 'w.'C  FI  GURU  A ,*r  PAST  PROGR  AN/RE?  ARA3LE  GH  i c.R  AT  I (Ni 
I DIDO  J ATA  RETRIEVAL  PROGRAM 

I 020  C 
I v)3-C 

104:W#S,R(SL)  *,8,13 

I .252$  s MSG  I :MUST  riE  RJN  AFTER  S.JJM-3  .VI  631'  PLEASE 
I 60 Si  IDENT*  -»P  I 19*, Ar ITSL  77  \ BRANTLEY  LORE  'AS 
I ’J  70S  * ( )PT I ON  : FOR  TP  A ii  , 0 " AP 
I 060$  : FORTY : NFOR ''  ,.iLNO 

I 090  0 

110  )0  PROGRAM  PAST  PROG RAM/ RE? ARABLE  GENERATION 

II  IOC 

i 1 200  THIS  PROGRAM  SEARCHES  THROUGH  EACH  5U A . TEE'S  GAT  A 
II 300  TO  RECORD  THE  PAST  OPERATIONAL  PROGRAM  \NJ 
11400  REPARABLE  GENERATIONS  FOR  EACH  MSN  IN  THE  DATA 
I 1 SCO  LIST. 

I I 600 

I 1 70  CHARACTER  F I IN*  I I ( 46)  , FN*  I I , DUN  HY*5(  24) 

I I 30  DIMENSION  PSTPEG( 46 , 35 ) , RPJN  ( 46 , 35  ) ,OFMDE  (4.5,  35  ) , 

I 190&FRCST(46f  35) 

1 2 50  DIMENSION  PSTPRG(  46, 35)  , RPGM  ( 46 , 35)  ,RGPP(2.4)  ,RG??U  -') 
1 2 1 PC 

12200  INITIALIZE  THE  TAPE  DEVICE  TO  ,T2  READ  FRO M IN  THE 
12300  DATA  BANT . 

I 24  JC 

12.50  INT=  I 2 
1 2600 

12/00  SET  THE  VALJH  FOR  THE  VARIABLE  REPRESENTING  THE 
12820  LAST  TAPE  DEVICE  IN  THE  DATA  SANK  TO  ;sE  READ  r-ROH 
12920  0,.  THIS  RUN. 

I 3000 

1315  IF  ( PSTPRG( 1 1 J) . 2Q. -I  . .OR . RPGH ( I , J ) ,E0. - 1 . ) 

I 32  0’,M5N(  I )='  )' 

133,5  IEN  = 1 2 
1 34  'C 

13500  READ  IN  THE  QUARTER  AND  NUMBER  OF  MSN. 

1 360  C 
1370  1=35 
1332  N J.'.‘=46 
1 390  0 

140 ’0  INITIALIZE  THE  PAST  PROGRAM  OR  SEPARABLE  GENEK'ATI 
14  100  TO  -1.  TO  PROVIDE  INDICATORS  FOR  JN  LATCHED  DSN'S. 

I 4 200 

143  : DO  I 230  N SN  = I ,4 6 

144  - DO  1 1 2.1  J=1  , M U M 

145  5 OFMDR  ( J oil , 36-1  ) =0 r ’ADR  ( ; i SN  ,23-1) 

I 4 , > ' RPG.'I  ( J,D  ) =-  I . 

1 4 7^0 

143  OF  HDRUSN,  ?2-I  ) = (0F  'DRC.ISM , 23- I ) +OFMDR  ( NSN , 2 1 - I ))/2. 
14900 

1520  OF’'  DR  ( N SN  , 7 ) =<  5 F PDR  ( .<  S. I , 5 ) = ( OF  R ( <3N,  3)  + 

I 5 I 'AOFMDR  ( NS!. , 5 ) )/2. 

1520  1 IP"  FOR  "AT  ( 14,  At  1 ) 

I 53  • OF  HOP  ( I3N,  5)  = 'F  'DR  (NSN,  I ) 


54  M 

READ  ( 1 >,  1 

-10, E 

i )=4 

j 3 

))PSLPRG 

( J 

,K) 

3 

RP 

DM  ( J , '< ) 

55" 

OFMOPCNSN, 

K)=. 9999 

560 

1110  CO:  IT  I 

NUE 

573 

FRO  ST  ( ..  SN  , 

J )=0F 

' >R  ( 

. i 5 

N,  J 

) +P  S'  7 

PR 

G ( . « 

5 

• ’ * 

J+l  ) 

580 

I =.  i SN 

590 

IF  ("SKI) 

.ED.' 

V ) 

GO 

TO 

1 2U" 

6 DOC 

6 1 3 

c 

62  3 

JO  1493  IN 

d:<=i  , 

260 330 

630 

GO  TO  114  0 

6 4 DC 

(j  50 

G DETER  MIN 

E IF 

THE 

LA 

ST 

AVAIL 

Ad 

LE 

T 

AP 

•*  1 * a or  . . ~ LT 

...  .MO  i > El  i_. 

663C 

6/3 

1135  INT=I 

NT+1 

680 

IF  ( INT.G'T 

. IEN) 

GO 

TO 

3 2 

m 

593 

1143  READ 

( INT, 

1 1 70 

» 

N D= 

1135) 

pT 

,0L 

A 

55 

,Fi.,(DUV.v 

700 

l)  , IMNTH=1  , 

24) 

710 

1170  FORMAT  (12 

0 V 

♦ t 

1 + 

, A 1 

i ,2X, 

24 

A 5 ) 

723  IF  (CLASS-6615)  1140,1171,1143 

730  1171  GO  TO  ( 1 1^3,  1 1 4 C , 1 20  3,  I 1 4(  ,1  14  0,  I 144,  1 1 4 A , I 1 4 J , 

74C.il  MS,  1 140,  1 200)  ,37 

750C 

760C  MATCH  THE  ITEM  READ  TO  THE  3ASIC  JATA  LIST 
77^C  PREVIOUSLY  SAVED.  Ir  TEE  STOCK  hUMsERS  MATO*, 

78 CC  INCREMENT  THE  COUNTER  (XNTR)  TO  TRACK  THE  NU'ViER 
79) C OF  DATA  POINTS.  IF  THE  STOCK  NUMBER  READ  DOESN'T 
6 0CC  MATCH  ANY  IN  THE  DATA  LIST,  READ  ANOTHER  j'U"BER. 

3 1 "C 

8 20  I 2 DC  DO  1473  J=l ,NUH 
8 30  IF  ( FN  .ME.  F I IN  ( J ) ) GO  TO  14  70 

840  kmtr=k:itr+i 


.'5  60  C 
370C 


SEPARATE  THE  MONTHLY  DATA  INTO  D'J AFTER 5 
THE  APPROPRIATE  MON  TIL Y JATA  TOGETHER. 


880  C 

890  .10= -2 


900  DO  1463  QTR=Q-4,Q 
910  3GPPQ  (OTR  >=3  . 

92'--  MQ=?«+3 

93D  DO  146J  ,MNTH=MQ,  "0  + 2 

940  IF  ( DU MMY ( MATH ) .EQ. •'***-’"> •')  GO  TO  1140 
950  DECODE  ( DU  "MY  ( MNTH)  , 1 453)  RGPP  ( ’/NTH) 

960  1+53  FORMAT (F 5.0) 

970  3GPPU (QTR ) =RGPPG ( QTR ) +3GPP( MMTH) 

930  146.0  CONTINUE 

993  IF  (RT.NE.03)  GO  TO  1 4 >5 


dY  ADDING 


2 P0D  3PGN ( J , OTR ) =RG PPO ( OTR ) 
2013  GO  TO  1463 


2.20  1455  PSTPRG ( J, QTR ) =3G PPO (DIR) 
2330  1463  CONTINUE 
2 04 0 IFCKNTR.LT. 2*NUM)  GO  TO  114 3 
2053  GO  TO  3000 
2 "60  14/3  CONTINUE 
2070  I 49M  CONTINUE 


4 6 


2 030  3000  00  3530  j = i,:jjm 
2090  PRINT  3010,  FSC( J > ,FI IN ( J) 

21  ) i 3 J I 3 FOP ’.‘AT  ( / , 1 4 , 1 X , A I I ) 

2110  (RITE  ( 1 7,332  4,EN0=4'O00)  r SC< J ) , F I IN < J) 

2123  3023  FOR WAT  ( 14 , A ! 1 ) 

2130  00  3200  0TR=1 , Q 

2140  PRINT  3100,  PSTPRGC J,QTR) ,RPGN< J,uTR) 

2150  3100  FORMAT  ( I I X,F7.  5,5X,F7.0> 

2160  TRITE  ( I 7, 31  20,  END=430O)  PSTPRC  ( J , J 2 R ) , R PGN ( J , 0 f R ) 
2170  3120  FORMAT  (F7.0,F7.0) 

2130  3200  CONTINUE 

2190  3530  CONTINUE 

2200  4000  COOT INUE 

2210  STOP 

2220  HMD 

2230S : EXECUTE 

224  3$ * LIMITS  *6  0,4  OX 

?25'-'$  s'fAPE  : I 2,  XI  0 , , 7 I 943  , , 77-4  00323 
226' $ s l A PE : : 6 , X2D , , 7 62  I 0 , , OATAF  I LE  5 
2270$  : FFILE*  1 7,NOSL2N 

223'M$:TAPEs  I /, AID, ,7361  3 , , OATAF  I LE9/P  I..3 
229'/$  jEMDJOO 


r 


(T00C  FIGURE  A.  5 DATA  PREPARATION  PROGRAM 

01  0C 

020C 


033 ##N , R( SL ) 

040$ : I DENT* WP 1193,  AFIT/SL  77 A E RANT LEY  LORE MAN 
353$  * OPTION *FORTR AN .NOMAP 
060$ » FORTY :NFORM,.iLMO 
373C 

0 SO C PROGRAM  DATA  PREPARATION 
390C 


1 00 C T^IS  PROGRAM  READS  THE  MASTER  STOCK  NUMBER,  OF" 

I13C  DEMAND  RATE,  BASE  REPARABLE  GENERATIONS,  AND  PAST 
I20C  PROGRAM  INFORMATION  FROM  T«E  FILE  AND  PERFORMS 
I 330'  CALCULATIONS  NECESSARY  TO  TRANSFORM  THEM  INTO  A 
I 43 C DATA  STREAM  FOP  ANALYSIS. 

1 50C 

163  CHARACTER  F I IN* 11 (46) , MSN  1 *1 5 (46) 

173  INTEGER  FSC(4->) 

133  DIMENSION  PSTPRGC 46 , 35 ) ,RPGN(46,35) ,OF’QR (46 , 35 ) , 

1 93AFRCSTC46, 35) 

2 33  C 

213C  READ  THE  MSN,  PAST  PROGRAM,  AND  REPARABLE 
223 C GENERATION  DATA  FROM  THE  DATA  TAPE. 

233C 

243  DO  1340  1=1,46 

250  READ  ( I 2, I O03,END=4  303)  F SC ( I ) , F I IN ( I ) 

263  I J00  FORMAT  (14, All) 

273  DO  1 030  J=l ,35 

283  READ  ( 1 2 , l 3 1 0 , EN Q=4Q 33 ) PSTPRC- ( I , J ) , RPGN  ( I , J ) 

290  1310  FORMAT  (2F7.0) 

300C 

310C  IF  PAST  PROGRAM  OR  REPARABLE  GENERATION  INFORMATION 
323 C IS  NOT  AVAILABLE  (EQUALS  -I.)  THEN  SET  THAT  FSC 
330 C EQUAL  TO  ZERO. 

343  C 

353  IF  ( PSTPRG( I,J) .EQ.-I . ,OR.RPGN( I , J) .EQ.-l .)FSC( I)»* 

363  1333  CONTINUE 
370  1343  CONTINUE 
330  DO  I 080  L = l ,46 
39  3 C 

4 03C  READ  THE  MSN  AMD  OF M DEMAND  DATA  FROM  THE  DATA  TAPE. 
4 I 0 C 

320  READ  ( 1 3, I353,END=4003)  MSNI(L) 

^ 33  1 053  FORMAT  (A15) 

443  DO  1073  K= I , P2 

450  READ  ( 1 3,  1060,E’ND=4  303)  OFMDP(L.K) 

Z63  1060  FORMAT  (F6.4) 

470  I 373  CONTINUE 


1480  1^30  CONTINUE 
1 490C 

1500C  FOLLOWING  SEGMENT  OF  PROGRAM  ‘5 
I 5 I 3 C OF  OF M DEMAND  DATA  NITH  THE  35 
I523C  PROGRAM  AND  REPARABLE  GENERATT 
I530C  MISSING  QUARTERS  OF  OFM  3EMAND 


SHES 


. J a 


DAF  \ 


THE  22  QUART 
TERS  OF  PAST 

ERE  ST'ULAT 


03 


48 


I 54  U C l3 Y SrRAIGHTL I •'£  AVERAGING  THE  VALUES  ON  EITHER  SITE. 
1 5500  FORECASTS  WERE  THEN  COMPUTED  3Y  MULTIPLYING  THE  Or" 
I560C  DEMAND  RATE  dY  THE  PAST  PROGRAM  Or'  THE  TEXT  QUARTER . 
I 5 7 0 C 

1560  DO  1280  MSN=I ,46 
1590  DO  1100  1=1,14 

1 6 00  OFMDRCiSM  , 36-  I )=()FMDR(N  SN  , 23-  I ) 
i6i0  Moo  a)NTi:.’UE 
1623  j:<=o 

I 630  JO  1 200  1 = 1,13,  2 
1640  J'<=JK+I 

1 650  OFMDP (NSN , 22-  I ) =(  OFMDR(NSN  , I 0-JK) +Or"DR  (.<  SN , 

1 66GA9- JK ) ) /2  . 

167  0 IiID=(  22-1  -5 ) /2 

1630  OFMDP ( N SN , 2 1 - I ) =0 FM  JR  (-1  SN,  I ND  ) 

1690  1200  CONTINUE 

1 700  OFMDP  (.1  SN , 7 ) = ( OFM DR  (N  SN  , 2 ) +OF"JR  (.■  3.1 , 1 ) )/2. 

1710  OFMDP ( MSN , 6 ) =OFMDR ( MS J,7) 

1 720  OFMDP (N SN , 5 ) =OFMDR (0  SN , 1 ) 

1730  DO  1235  *f=l,4 
1 740  OFMDR  (MSN,  !<)=.  9999 
1750  1235  CONTINUE 
1760  DO  1210  J=5,34 

1 770  FRCST(NSN,J)=OFMDR(NSN, J)*PSTPPG(NSN, J+l  ) 

I780  1210  CONTINUE 
1790  I=N  SN 
1 8 00 C 

18I0C  PRINT  AND  WRITE  TO  A PERMANENT  FILE  THE  OFM 
1 820C  DEMAND  RATE,  PAST  PROGRAM,  FORECAST,  AMD  SEPARABLE 
I 8 33  C GENERATIONS  FOR  EACH  QUARTER. 

1 640 C IF  THE  FSC  EQJALS  0 (NO  VALID  DATA),  TH  *N  SHI?  IT. 

1 850C 

1863  IF  (FSC(  D.53.3)  GOTO  123J 
1870  PRINT  1230,  FSCC I ) , F I IN ( I ) 

1 330  1 230  FORMAT  (//, 1 3X,3HMSN,2X, 14, A 11,/) 

1 a93  WRITE  (17.124U)  FSC(  I ) , FI  IN  ( I ) 

1 9 30  1 24  3 FORMAT  ( 14,  A!  I ) 

1910  PRINT  1245 

1920  1245  FORMAT  ( 24  X,  4 HP  ASF,  2 3X  , 9HRERAR  A^LE , / , I X , 

1 9 30  5<  7 '-/QUARTER,  5X , 5H0FMJR , 5X  , 7 HP  POOR  AM,  5X,  3HFCRSCASF  , 
I940A5X, 1 1HGEHERATI0NS,//) 

1 950  DO  1270  J=1 ,35 

1960  PRINT  I 250,  J,OFMDR( I , J) ,PSTPRG( I , J ) ,F  RCST( I , J ) , 

1 9701.RPGN  ( I,  J) 

1980  1250  FORMAT  (4X,I2,6X,F6.4,5X,F6.  ’ , 5X , F6 .4 , 3X ,F 5. 0 ) 
1990  NR  ITS  (17,1280)  J , OFMJR ( I , J ), PST?RG( I , J ), ERCSF ( I , J ) , 

2 O00&RPGN ( I , J ) 

2010  1263  FORMAT  ( I 2,  2X,  F6. 4 , 2X,  F6 .0, 2X,r  3 ,4 , 2X,  F6.  ) 

"023  1270  CONTINUE 
2030  1230  CONTINUE 
2.040  4 000  ONTINUE 
2C5r/  STOP 
2060  END 
2r  7'’Ss  EXECUTE 
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TABLE  A . 1 


D041  DATA 


Quarter 

OFM 

Demand  Rate 

Past 

Program 

Reparable 

Generations 

1. 

69-2 

★ 

1016 

28 

2. 

69-3 

★ 

936 

25 

3. 

69-4 

* 

914 

23 

4. 

70-1 

★ 

861 

30 

5. 

70-2 

. 0336 

816 

24 

6. 

7 0-3 

★ 

791 

21 

7. 

7 0-4 

★ 

693 

32 

8. 

71-1 

. 0314 

571 

24 

9. 

71-2 

★ 

580 

22 

10. 

71-3 

. 0357 

621 

33 

11. 

71-4 

* 

548 

21 

12. 

72-1 

.0380 

479 

17 

13. 

72-2 

★ 

497 

32 

14. 

72-3 

. 0433 

533 

15 

15. 

72-4 

* 

502 

17 

16. 

7 3-1 

. 0424 

493 

24 

17. 

73-2 

* 

577 

17 

18. 

73-3 

.0367 

809 

16 

19. 

73-4 

* 

928 

27 

20. 

74-1 

. 0301 

981 

11 

21. 

74-2 

★ 

663 

14 

22. 

74-3 

. 0243 

652 

10 

23. 

74-4 

. 0235 

435 

11 

24. 

75-1 

.0218 

285 

10 

25. 

75-2 

. 0213 

364 

10 

26. 

7 5-3 

. 0222 

383 

11 

27. 

75-4 

. 0202 

392 

7 

28 

76-1 

.0247 

384 

15 

29. 

76-2 

. 0266 

342 

12 

30. 

76-3 

. 0283 

379 

8 

31. 

76-4 

.0293 

384 

15 

32. 

77-1 

. 0300 

342 

12 

33. 

77-2 

. 0339 

379 

8 

34. 

77-3 

.0335 

335 

11 

35. 

77-4 

.0357 

324 

13 

*Indicates  missing  data. 
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TABLE  A. 2 


DATA  AFTER  PREPARATION 
MSN  6615000735893 


PAST 

REPARABLE 

iRTER 

0FMDR 

PROGRAM 

FORECAST 

generations 

1 

0,9999 

1016. 

0 . 

28  . 

2 

0 . 9999 

9 36  . 

0 . 

25  . 

3 

0.9999 

914. 

0 . 

23. 

4 

0.9999 

861  . 

0 . 

30  . 

5 

0.0336 

816. 

26.5776 

24  . 

6 

0 . 0 325 

791. 

22.5225 

21  . 

7 

0 . 0 325 

6 93  . 

18.5575 

32. 

8 

0.0314 

571. 

18.2120 

24  . 

9 

0 . 0 336 

5 8 0 , 

20.8346 

22  . 

1 0 

0 . 0 35  7 

6?1  . 

19.5636 

33  . 

1 1 

0 . 0 369 

5 4 8 . 

17.6512 

21  . 

12 

0 . 0 38  0 

4 79  . 

18.8860 

17. 

13 

0 . 0 4 0 7 

4 97  . 

21.6665 

32  . 

1 4 

0.0433 

5 33  , 

21.7366 

15  . 

15 

0 . 0 42  9 

5 0 2 . 

21.1251 

17. 

16 

0.0424 

4 93  , 

24.4648 

24  . 

17 

0.0396 

577  . 

31.9960 

17  . 

1 8 

0.0  367 

8 0 9 , 

34.0576 

16  . 

19 

0 . 033  4 

92  8 . 

32.7654 

27  . 

20 

0.0301 

981  . 

19.9563 

11  . 

21 

0.0272 

6 6 3 , 

17.7344 

14  . 

22 

0.0  24  3 

652  , 

10.5705 

10. 

23 

0 . 0 235 

4 35  . 

6.6975 

11  . 

24 

0.0218 

2 85  . 

7.9352 

10  . 

25 

0.0213 

36  4 , 

8.1579 

10  . 

26 

0 , 0 222 

3 83  , 

8.7024 

11 . 

27 

0 . 0 2 0 2 

3 92  . 

7.7568 

7 . 

28 

0 . 0 2 4 7 

3 84  . 

8,4474 

15  . 

29 

0 . 0 26  6 

3 4 2 . 

10.0814 

12  . 

30 

0.0  283 

3 79  . 

1 0 . 8 672 

8 . 

31 

0.0293 

3 8 4 . 

1 0 . 0 2 0 6 

15. 

32 

0.0300 

3 4 2 . 

11.3700 

12  . 

3 3 

0,0339 

379  . 

11.3565 

8 . 

34 

0 . 0 335 

335  . 

1 0 . 8 54  0 

11  . 

35 

0.0357 

3 2 4 . 

0 . 

13  . 
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FOR  HO AST 


)34  I 


1 030 0 THIS  PROGRAM  ANALYZE 
I ..'40C  T"E  ERROR  ASSOCIATED 

i os.' c lead  times  one,  ;..o, 

I FOv'C 

I 079  CHARACTER  MSN*  I 3 
I DIMENSION  FRCoIK  35)  ,Rpi 


. 90C 


— pi.  r ' 


Of)  ' 


T ... 


THE  L>04 1 u’A:A  A..J  Of) 
i IT':  :£ ACM  FORECAST  TF 
T ’ ' , A.  i . ' r./  J P . 


N(  35  ) , EPR.;04  I (40,.. ) 


li  ' : C REA.)  FROM  A FILE  THE  MASTER  STOCK  lU.’lSIR, 
I I I PC  REPARABLE  GENERATIONS  FOR  EACH  OUATIER. 
i I 2 ’ C 

1130  C \LL  ATTACK  I 5,  "77  V34/ )AT)"4  1 ',1  ,,) 

1140  READ  (16,900)  LN , "SN 
I I 50  9 -X1  FOR  MAT  (V) 

I 1 60  DO  1 000  1 = 1 , 35 

I I 70  READ  (16,900)  LN, FRCST ( I ) , RPGN( I) 

1130  13  00  CONTINUE 
1 1 90  DO  I I on  1=3,34 
1 300  DO  1090  J=1  ,4 
13  10  i;iD  = I+J 

1320  IF  (IND.GT.35)  'JO  TO  I .’SO 


23.’ 

04  1 (I , ,J  ) = 

FROST 

( 

I) -PPG. 

.i(  I 

NO) 

34.) 

(j(;  77)  j ^91 

25D 

1 330 

CONTINUE 

262 

1 J9.) 

ONTIN'JE 

270 

1 1 NO 

CONTINUE 

2 6 '2' 

pr  1.: 

TIM) 

290 

Ml) 

FORMAT  ( 

34  X , 11 

ABLE  A 

< II 

, // ,24  X , 

300’ 

O' COM 

PUT  ED  D34 

1 FOR 

CASr  E 

I »/////) 

310 

PR  IN' 

T 112) 

320 

1 1 23 

FORMAT  C 

21  X, 4 

( 

2X , 3HL 

had 

TIME)  ) 

330 

PR  IN 

T 1130 

340 

1 1 3 J 

FORMAT  ( 

27  X,  1 

1 ,9X,  i 

M2, 

9X,  I'M, OX, 

350 

PRIi." 

T 1 1 35 

360 

1 1 35 

FORMAT  ( 

1 5 X , 7 

!■ 

j U A 7 i .1 

R , 3 

X , 5 HERR  0 ? , 

372 

DO  i 

1 60  1=5,3 

4 

330 

PRIM 

— 1 I 4 • 3 T 

ERROR 

i c 1, 1 ) 

, E RDM  '41  (1,2) 

390- 

IERRQ 

04  1 C I , 3)  , 

ERR30 

4 

1(1,4) 

4 JO 

1140 

FOR’; AT  ( 

1 SX,  I 

0 

,4 (4X, 

r j • 

1 ) ,/) 

4 | 0 

1 1 50 

CONTINUE 

4 20 

1 1 60 

CONTINUE 

4 30 

STOP 

44-  END 


1 F4 ) 

3 ( 5 X , 5 


:P  J D ;S 
OJ  J'f 


IS  )Ai 


. )R ) , //> 
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APPENDIX  B 

TIME  SERIES  FORECASTING  DATA  PREPARATION 


w 


The  forecasting  method  used  was  that  of  the  Honey- 
well Corporation  TCAST  computer  program  in  the  AFIT  CREATE 
Library.  This  program  reviews  the  exhibited  characteris- 
tics of  a variable  and  describes  its  behavior  by  perform- 
ing an  intrinsic  analysis  on  discrete  data.  In  the  analysis 
each  data  point  is  split  into  base  series,  trends,  seasonal, 
cyclical  and  erratic  components.  Base  series  are  known 
historical  data  values  which  are  considered  constant  over 
given  periods  of  time;  trend  refers  to  the  long-term  growth 
or  decay  of  the  system;  seasonal  refers  to  regular  annually 
recurring  periodic  fluctuations  which  are  predictable; 
cyclical  is  a wave- like  fluctuation  about  a trend  which  is 
not  constant  nor  predictable;  erratic  refers  to  the  unex- 
plained or  random  variation  that  occurs  for  unknown  reasons. 
After  the  analysis  on  each  data  component  the  results  are 
combined  and  the  forecasting  model  developed  generates  a 
forecast . 

The  forecast  generation  followed  five  objectives: 

1.  Evaluated  trends  in  the  historical  data. 

2.  Dampened  the  erratic  variations. 

3.  Responded  to  true  changes  in  the  physical 
process  as  they  occurred. 

4.  Projected  trends  in  future  time  periods. 
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5.  Used  program  generated  optimum  forecast  param- 
eters to  accomplish  the  previous  four  objectives. 

A "signal,"  the  deterministic  part  of  the  time 
series,  consists  of  the  long-term  trend,  cyclic  and  seasonal 
effects  and  this  signal  is  separate  from  the  random  vari- 
ation. This  signal  and  the  random  variation  projections 
are  used  together  to  determine  a future  process  value. 
Logically  these  steps  were  followed  to  develop  a signal: 

1.  Display  the  data  and  visually  determine  any 
obvious  characteristics. 

2.  Remove  any  apparent  growth  or  decay  trends. 

3.  Remove  any  significant  cyclical  components. 

4.  Evaluate  the  residual  trends  of  the  original 
corrected  data  and  determine  their  type  and  order  of 
growth  or  decay. 

5.  Fit  the  residual  trends  to  an  appropriate 
program  method:  exponential  smoothing. 

6.  Forecast  for  the  future  period  using  the 
appropriate  programmed  time  series  forecast  model. 

7.  Continue  to  update  the  model  as  new  data 
becomes  available. 

As  updating  continues  various  alphas  are  selected  within  the 
TCAST  Program  to  give  the  smallest  mean  absolute  deviation 
(MAD)  of  the  forecast  errors  for  each  smoothing  type. 
Generally,  alphas  will  be  between  .01  and  .3  for  first 
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order  smoothing  and  . 3 to  one  for  second  or  third  order 
smoothing.  After  smoothing  has  been  accomplished  the 
residual  left  is  due  to  random  variations  and  these  would 
produce  errors  in  any  forecast  model. 

Forecasts  for  this  study  were  generated  by  trun- 
cating the  TCAST  program  through  the  use  of  ASK  function. 
Accordingly,  all  extraneous  data  was  sent  to  a file  dump 
while  only  the  forecast  information  was  asked  for  on  the 
teletype,  thereby  decreasing  the  actual  time  for  the  te_e- 
type  operation  of  TCAST.  The  actual  operation  of  TCAST 
was  easily  accomplished  by  following  the  instructions  that 
can  be  obtained  from  AFIT/CREATE  Library. 

The  forecasts  generated  were  obtained  by  using  the 
first  eight  data  points  to  establish  a computational  base 
for  the  first  forecast,  with  additional  data  items  input 
to  continue  forecasting  into  future  periods.  Forecasts 
were  generated  for  lead  times  (number  of  periods  in  the 
future  being  forecasted)  of  one,  two,  three  and  four. 
Specifically,  for  lead  time  one,  forecasts  began  in  the 
ninth  quarter  and  continued  to  the  35th  quarter;  for  lead 
time  two,  forecasts  began  in  the  tenth  quarter  and  con- 
tinued to  the  thirty-fifth;  for  lead  time  three;  fore- 
casts began  in  the  eleventh  quarter  and  continued  to  the 
thirty- fifth  quarter;  and  for  lead  time  four  fore-casts 
began  in  the  twelftn  quarter  and  continued  to  the  thirty- 
fifth  quarter.  One  complete  TCAST  run  had  to  be  made  for 
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each  needed  forecast;  i.e.,  the  first  run  included  the 
first  eight  data  points  to  form  the  data  base  to  forecast 
for  the  ninth  period,  the  second  run  included  the  data 
base  plus  the  influence  of  the  ninth  periods  actual  usage 
to  forecast  for  the  tenth  period,  and  so  forth.  Each  run 
continued  to  regenerate  forecast  for  the  past  periods  for 
which  a forecast  had  already  been  extracted  but  these  were 
disregarded  and  only  the  last  period  forecasted  on  each 
run  was  extracted. 

After  the  forecasts  were  obtained  for  each  quarter 
the  type  smoothing  (first  order,  second  order,  or  third 
order)  was  recorded,  as  was  the  alpha  (optimal  weighted 
value) . Following  this  the  actual  reparable  generations 
were  subtracted  from  the  forecasted  values  and  this  differ- 
ence represented  the  error  for  each  quarter.  The  error 
was  rounded  to  one  decimal  place  and  both  positive  and 
negative  values  were  obtained.  (See  Figure  B.l  through 
B.4.)  These  differences  represented  the  error  distribu- 


tions data  base. 


TABLE  B.l 


TCAST  DATA  FOR  LEAD  TIME  ONE 


Qtr. 

Order 

Smoothing 

Smoothing 

Constant 

Reparable 

Forecast  Generations 

Forecast 

Error 

8 

3 

.2 

25.453 

24 

9 

2 

.8 

19.429 

22 

3.5 

10 

1 

. 3 

27.253 

33 

-13.6 

11 

3 

.1 

24.953 

21 

6.3 

12 

3 

.1 

22.370 

17 

8.0 

13 

1 

.1 

26.092 

32 

-9.6 

14 

3 

.1 

21.706 

15 

11.1 

15 

3 

.1 

19.830 

17 

4.7 

16 

1 

.2 

22.536 

24 

-4.2 

17 

3 

.1 

18.911 

17 

5.5 

18 

3 

. 1 

17.407 

16 

2.9 

19 

3 

.1 

19.550 

27 

-9.6 

20 

3 

. 1 

16. 516 

11 

8.6 

21 

3 

. 1 

15. 026 

14 

2.5 

22 

3 

. 2 

10.008 

10 

5.0 

23 

3 

.2 

8.8349 

11 

-1.  0 

24 

3 

. 2 

7.8006 

10 

-1.2 

25 

2 

. 2 

8. 6828 

10 

-2.2 

26 

2 

.2 

8.6299 

11 

-2.3 

27 

3 

.1 

6.7158 

7 

1.6 

28 

2 

. 2 

9.3022 

15 

-8.3 

29 

1 

.6 

12.156 

12 

-2.7 

30 

3 

.1 

7.2936 

8 

4.2 

31 

2 

. 2 

10. 514 

15 

-7.7 

32 

2 

. 3 

12.162 

12 

-1.5 

33 

2 

. 2 

9.3827 

8 

4.2 

34 

1 

.2 

7.7847 

11 

-1.6 

35 

15 

-5.2 

TABLE  B . 2 


Qtr . 

Order 

Smoothing 

Smoothing 

Constant 

Forecast 

Reparable 

Generations 

Forecast 

Error 

9 

1 

. 3 

25. 985 

10 

3 

.1 

24.098 

33 

-7.0 

11 

3 

. 1 

26. 659 

21 

3.  9 

12 

3 

.1 

24.756 

17 

9.7 

13 

3 

.1 

21.929 

32 

-7.3 

14 

3 

.1 

24.652 

15 

6.9 

15 

3 

.1 

21. 242 

17 

7.7 

16 

3 

. 1 

19.211 

24 

2.  8 

17 

3 

. 1 

19.950 

17 

2.2 

18 

3 

.1 

18.281 

16 

4.0 

19 

3 

.1 

16.672 

27 

-8.7 

20 

3 

. 1 

19.081 

11 

5.  7 

21 

3 

. 1 

15. 781 

14 

5.1 

22 

3 

.1 

14.195 

10 

5.8 

23 

3 

.1 

11.684 

11 

3.2 

24 

3 

. 1 

10. 096 

10 

1.7 

25 

3 

. 1 

8. 5976 

10 

. 1 

26 

3 

. 1 

7. 5072 

11 

-2.4 

27 

3 

.1 

7.0772 

7 

. 5 

28 

3 

. 1 

5. 5248 

15 

-7.9 

29 

3 

.1 

7. 0390 

12 

-6.5 

30 

2 

. 2 

9.2174 

8 

-1.  0 

31 

2 

. 2 

7. 9214 

15 

-5.8 

32 

2 

.2 

10.175 

12 

-4.1 

33 

2 

.2 

10.491 

8 

2.2 

34 

1 

. 4 

9.6932 

11 

- .5 

35 

13 

-3.3 
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TABLE  B . 3 


TCAST  DATA  FOR  LEAD  TIME  THREE 


Qtr . 

Order 

Smoothing 

Smoothing 

Constant 

Reparable 

Forecast  Generations 

Forecast 

Error 

10 

2 

. 5 

25.117 

11 

1 

. 8 

22. 638 

21 

4.1 

12 

2 

. 2 

27.594 

17 

5.6 

13 

3 

.2 

23.838 

32 

-4.4 

14 

3 

.2 

17.109 

15 

8.8 

15 

3 

. 2 

26.159 

17 

. 1 

16 

3 

. 2 

17.179 

24 

2.2 

17 

3 

. 2 

13.819 

17 

.2 

18 

3 

.2 

18.207 

16 

-2.2 

19 

3 

. 2 

15.177 

27 

-2.2 

20 

3 

.2 

12.870 

11 

4.2 

21 

3 

. 2 

21.246 

14 

-1.1 

22 

3 

. 2 

12.779 

10 

11.2 

23 

3 

.2 

10. 619 

11 

1.8 

24 

3 

.2 

6.3872 

10 

. 6 

25 

3 

. 2 

5.2920 

10 

-3.6 

26 

3 

.2 

4.3943 

11 

-5.  7 

27 

3 

. 2 

4.4559 

7 

-2.6 

28 

3 

.2 

5.8945 

15 

-10.  5 

29 

2 

. 2 

5.1863 

12 

-6.1 

30 

2 

. 2 

8.0304 

8 

-2.8 

31 

2 

.2 

8.6894 

15 

-7.0 

32 

1 

. 8 

8.8621 

12 

-3.3 

33 

1 

. 8 

13.772 

8 

. 9 

34 

1 

. 4 

10.338 

11 

2.8 

35 

13 

-2.7 
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TABLE  B . 4 


TCAST  DATA  FOR  LEAD  TIME  FOUR 


Order 

Smoothing 

Reparable 

Forecast 

Qtr.  Smoothing 

Constant  Forecast 

Generations 

Error 

11 

3 

. 3 

8. 0450 

12 

3 

.3 

18.625 

17 

-9.  0 

13 

3 

. 2 

30.163 

32 

-13.4 

14 

3 

. 2 

23.476 

15 

15.2 

15 

3 

. 2 

15.709 

17 

6.5 

16 

3 

. 2 

26.359 

24 

-8.3 

17 

3 

. 2 

15.927 

17 

9.  4 

18 

3 

. 2 

12.164 

16 

- .1 

19 

3 

.2 

17.436 

27 

-14.8 

20 

3 

. 1 

16.967 

11 

6.4 

21 

3 

. 2 

11.475 

14 

3.0 

22 

3 

. 1 

18.112 

10 

1.5 

23 

3 

. 2 

11.548 

11 

7.1 

24 

3 

. 2 

9.1894 

10 

1.  5 

25 

3 

.2 

4.4506 

10 

- .8 

26 

3 

. 2 

3.4062 

11 

-6.6 

27 

3 

. 2 

2.5908 

7 

-3.6 

28 

3 

. 2 

2.8846 

15 

-12.4 

29 

3 

. 2 

4.7567 

12 

-9.1 

30 

3 

.2 

2.5082 

8 

-3.2 

31 

3 

.1 

5.0361 

15 

-12.  5 

32 

3 

.1 

5.5723 

12 

-7.0 

33 

2 

. 2 

6.7258 

8 

-2.4 

34 

1 

. 4 

8.5262 

11 

-4.3 

35 

13 

-4.5 
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APPENDIJ 
STATISTICS  COI- 


The  statistics  for  the  error  distributions  were 


computed  using  the  LSS/SIMFIT  program.  The  program  fitted 
the  error  data  points  to  an  internally  determined  distribu- 
tion type,  and  then  provided  the  type  distribution  derived 
and  the  statistics  as  shown  in  Figures  C.l  through  C.8. 
These  statistics  were  then  used  to  perform  Student's  t 
tests  on  the  null  hypotheses. 

The  SIMFIT  Program  gave  the  uniform  distribution  as 
the  predominate  distribution  fitting  the  data  (smallest 
maximum  error)  for  five  of  the  eight  error  distributions. 
However,  the  histogram  provided  by  SIMFIT  empirically  dis- 
played the  data  in  symetrical  form.  The  normal  distribu- 
tion's maximum  error  was  larger  than  the  Uniform's,  how- 
ever, in  testing  for  the  Uniform  the  larger  errors  were  not 
sufficient  to  reject  the  Uniform  distribution  at  the  90 
percent  confidence  level.  But  since  empirical  evidence 
must  be  weighted  more  heavily  than  theoretical,  the  normal 
distribution  was  accepted  as  being  the  predominant  distri- 
bution, which  supported  the  previous  assumption  that  the 
forecast  error  was  normally  distributed.  This  can  be  more 
clearly  seen  in  the  graphs  in  Figures  C.9  through  C.16. 

' The  slightly  larger  variance  of  the  D041  forecasts 
can  be  easily  seen  in  the  graphs  (Figure  C.17  through  C.24) 
which  plot  the  forecast  error  over  time. 
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THIS  IS  THE  ERROR  DISTRIBUTION 


SAMPLE  STATISTICS 

NUMBER  OP  VALUES  IN  LIMITS  = 30.000 

MINIMUM  VALUE  OF  ALL  DATA  = -13.100 

MAXIMUM  VALUE  or  ALL  DATA  = 21,800 

MODE  Or  DATA  = -3.750 

ARITHMETIC  MEAN  = 0.517 

S T f)  . DEVIATION  = 7.476 

VARIANCE  = 55.897 

VARIANCE  / mean  = 108.188 

ESTIMATED  SAMPLE  MOMENTS  JPOUT  THE  MEAN: 

3 276.4 

4 0.1167E  05 


A3 

0.696 

I r 

A 3 

IS  EQUAL  TO 

ZERO.  THE  DISTRIBUTION  IS  SYMMETRICAL 

A 4 

0.997 

I r 

A 1 

IS  NEGATIVE 

DISTRIBUTION 

Has  A FLAT  TOP 

i r 

A 4 

IS  Z F R O . DISTRIBUTION  is 

ALMOST  NORMAL 

I r 

A 4 

IS  POSITIVE, 

DISTRIBUTION 

has  a high  p r a ¥ . 

i r 

A NY 

PROBABILITY 

FUNCTIONS  OO 

NOT  FIT  DATA 

SUCRE  ST  THAT  YOU  TRY  THE  FOLLOWING  VALUES 
AMEN  = -13.1000 

U I 0 T H = 3.8778 

N C P !.  L S = 9 


Figure  C.l.  Statistics  Computations  for  D041 
Forecasts  in  Lead  Time  One 
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THIS  is  T M F 6;' POP  DISTRIBUTION 


S A M P l E STATISTICS 


NUMRI:D  OF  V4LIJFS 

I M LIMITS  = 

? 9 

.000 

MINI  MUM  ' A L ME  OF 

all  data  = 

- 1 4 

.800 

MAXIMUM  VALUE  OF 

AIL  DATA  = 

2 3 

.10  0 

MODE  OF  RAT*.  = 

0.750 

4P1TiImET!C  mean  = 

0 , B B 3 

S T l'  . DEVIATION  = 

7.908 

VARIANCE  = 

6?. 536 

VARIANCE  / “FAN  = 

70.84? 

ESTIMATED  SAMPLE 

moments  ABOUT 

THE 

MF  A N 

3 3 35  .3 

4 0.165  -IE  05 


A3  0.715 

IE  43  IS  EQUAL  TO  ZERO.  THF  DISTRIBUTION  IS  SYMMETRICAL 
A 4 1.536 

IF  44  IS  NFGAT1VE  DISTRIBUTION  HAS  A FLAT  TOP 

IE  4 4 IS  ZERO,  DISTRIBUTION  IS  4 1.  MOST  NORMAL 

IF  44  IS  POSITIVE,  DISTRIBUTION  HAS  A HIGH  PEAK. 

IF  ANY  PROBABILITY  FUNCTIONS  DO  NOT  TIT  DATA 
SUGGEST  THAT  YOU  TPY  TOP  FOLLOWING  VALUES 
APF‘i  = -14. 8000 

WIDTH  = 4.3563 

N C F L L S = 9 

Figure  C.2.  Statistics  Computations  for  D041 
Forecasts  in  Lead  Time  Two 
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THIS  IS  THE  ERROR  DISTRIBUTION 
SAMPLE  STATISTICS 

NIJMRER  of  VALUES  IN  LIMITS  = 28.000 
MINIMUM  VALUE  OF  ALL  DATA  = -14.400 
MAXIMUM  VALUE  OF  ALL  DATA  = 22.800 


MODE  OF  DATA  = 0.750 
ARITHMETIC  MEAN  = 1.650 
STD.  DEVIATION  = 8.662 
VARIANCE  = 75.031 
VARIANCE  / MEAN  i 45.473 


ESTIMATE^  SAMPLE  MOMENTS  AROUT  the  mean: 

3 534.1 

4 0.2II67E  05 

A3  0.868 

IF  A3  IS  EUUAL  TO  ZERO,  THE  DISTRIBUTION  IS  SYMMETRICA! 
A 4 0 . 9 4 9 

IF  A 4 IS  NEGATIVE  DISTRIBUTION  HAS  A FLAT  TOP 

IF  A 4 IS  ZERO.  DISTRIBUTION  IS  ALMOST  NORMAL 

IF  A 4 IS  POSITIVE,  DISTRIBUTION  HAS  A HIGH  PEAK. 

IF  ANY  PROBABILITY  FUNCTIONS  DO  NOT  TIT  DATA 
SUGGEST  THAT  YOU  TRY  THE  FOLLOWING  VALUES 
AMEN  i -14.4000 

WIDTH  = 4.4286 

NOFILS=  8 

Figure  C.3.  Statistics  Computations  for  D041 
Forecasts  in  Lead  Time  Three 


THIS  IS  TUP  ERROR  DISTRIBUTION 


gamble  statistics; 


NUMBER  OF  VALUFS 

IN  LIMITS  = 

2 7 

.000 

MINIMUM 

V A 1 UE  OF 

ALL  DATA  = 

-11 

.200 

MAXIMUM 

VALUE  OF 

ALL  O A T A = 

2 4 

.10  0 

MOPE  OF 

) A r A = 

0.750 

ARITHMETIC  MEAN  = 

2.170 

STD.  DEVIATION  * 

0.069 

variance 

s 

8 2.740 

V A R I A N C E 

/ MEAN  = 

37  . 76  7 

E S T I 1 A TE 

D SAMPLE 

MOMENTS  APOUT 

THE 

MEAN 

6 14.0 

0 . 2 0 2 6 E 05 


AT  0.071 

[F  4.5  IS  EOUAL  TH  7ERO,  T H F DISTRIRUTION  IS  SYMMfTRICAl 
44  0.230 

IF  4 4 IS  NEGATIVE  DISTRIBUTION  HAS  A FLAT  TOP 

IF  \ A IS  2 E D 0 < DISTRIBUTION  IS  ALMOST  NORMAL 

IF  A 4 IS  POSITIVE.  DISTRIBUTION  HAS  A HIGH  PEAK. 

IF  ANY  ° R OH  A ij  1 1.  I T Y FUNCTIONS  DO  NOT  FIT  DATA 
SUGGEST  THAT  YOU  TRY  THE  FOLLOWING  VALUES 
AME'I  = -11,2000 

WIDTH  = 4.3500 

NCEI LS=  S 


Figure  C.4.  Statistics  Computations  for  D041 
Forecasts  in  Lead  Time  Four 
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THI‘>  IS  T II F ERROR  DISTRIBUTION 
SAMP!  E STATISTICS 


NUMBER 

or  VALUES 

I H 

L I M I rs  = 

27 .0  0 n 

MINIMUM 

V A LUC 

OF 

Al  L 

DATA  = 

-13.600 

MAXIMUM 

V A 1.  u E 

OF 

ALL 

DATA  = 

11.10  0 

mode  or 

DATA 

= 

S . ? 5 0 

APIT'IMCTIC  “<F  AN  = -0*096 

STD.  DEVIATION  : 6 . ?31 

VARIANCE  s j 0 . 0 ? 7 

VARIANCE  / MEAN  = -403.  ? 0 7 

ESTI^ATFO  SAMPLF  MOMENTS  APOUT  THE  MEAN: 

3 -70.11 

4 3 7 9 9, 


A3  -0.307 

IE  AS  IS  EQUAL  TO  ZERO.  T H F DISTRIBUTION  IS  SYMMETRICAL 


A 4 

-0 

.640 

I F 

A 4 

I S 

NEGATIVE 

DISTRIBUTION 

HAS  A 

FLAT  TOP 

I E 

A 4 

I s 

7FRO,  DISTRIBUTION  IS 

ALMOST 

NORMAL 

I F 

A 4 

I s 

POSITIVE. 

D I STB  I BUT  I ON 

HAS  A 

HIGH  PEAK. 

I F 

AMY 

p 

ROHABILITY 

functions  DO 

NOT  FIT  DATA 

SIIO  or  ST  THAT  YOU  TRY  THF  FOLLOWING  VALUES 
AMF'I  - -13.6000 

WIDTH  = 3.0494 

NCELL5*  0 

Figure  C.5.  Statistics  Computations  for  TCAST 
Forecasts  in  Lead  Time  one 
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THIS  IS  Tii-  ERROR  DISTRIBUTION 
SAMPLE  STATISTICS 


number 

OF  VALUES 

If;  LIMITS  = 

26.000 

MINIMUM 

7 A L HE 

OF 

all  data  = 

-«  . 7 0 0 

MAXIMUM 

VALUE 

or 

all  data  = 

9.700 

MODE  OF 

DATA 

= 

0.750 

ARITHMETIC  MEAN  s D.269 

ST'!.  DEVIATION  = 5 . 3 0 B 

VARIANCE  = 2 S . 1 7 4 

VARIANCE  / MEAN  = 10  1.  A 4 7 

ESTIMATED  SAMPLE  MOMENTS  ABOUT  THE  MEAN: 

3 -21.13 

1 1115. 

A .3  -0.171 

IF  A3  IS  EQUAL  TO  2ERO,  THE  DISTRIBUTION  IS  STMMETRICA! 
A 4 -1.072 

IF  A 1 IS  NEGATIVE  DISTRIBUTION  HAS  A FLAT  TOP 

IF  A A IS  ZERO,  DISTRIBUTION  IS  ALMOST  NORMAL 

IF  AT  IS  POSITIVE,  DISTRIBUTION  HAS  A HIGH  PEAK. 

IF  AMY  PROBABILITY  FUNCTIONS  DO  NOT  FIT  DATA 
SUGGEST  THAT  YOU  TRY  T H R FOLLOWING  VALUES 
a m e ' = -fl,7onn 

W I n th  = 2.3590 

N C E L L S = 0 

Figure  C.6.  Statistics  Computation-  for  TCAST 
Forecasts  in  Lead  Time  Two 
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THIS  IS  TUT  ERROR  DISTRIBUTION 


SAMPLE  STATISTICS 


N U M 0 F P nr  VALUES 

IN  LIMITS  = 

25 

.000 

minimum  value  nr 

ALL  DATA  = 

-1  0 

.500 

MAXIMUM  VALUE  or 

ALL  DATA  = 

1 1 

.200 

modf  or  data  = 

0.750 

ARITHMETIC  MEAN  = 

-0.816 

STP.  DEVIATION  = 

5.231 

VARIANCE  = 

2 7.363 

VARIANCE  / M F A N r 

-53.029 

ESTIMATED  SAMPLE 

MOMENTS  AROUT 

THE 

mean 

3 

26.82 

4 

1 9 0 3 . 

A3 

0 . 

19  9 

I F 

A 3 

I s 

EQUAL  TO  ZERO.  THE  DISTRIBUTION  IS  SYMMETRICAI 

A 4 

-0  . 

2 4 2 

I F 

A t 

I s 

NEGATIVE  DISTRIBUTION 

HAS  A FLAT  TOP 

i r 

A 4 

I s 

ZERO.  DISTRIBUTION  IS 

almost  normal 

I r 

A 4 

I s 

POSITIVE,  distribution 

HAS  A HIGH  PEAK. 

i r 

ANY 

probability  functions  do 

NOT  FIT  DATA 

SUGGEST  THAT  YOU  TRY  THE  FOLLOWING  VALUES 

AMES  = -in. 5000 

WIDTH  = 2.8933 

NCE  I I.  S=  8 


Figure  C.7.  Statistics  Computations  for  TCAST 
Forecasts  in  Lead  Time  Three 
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THIS  15  THF  FRROR  DISTRlRIlTION 


SAMPLE  STATISTICS 

I 

I 


1!  !I»RTR  OF  VALUES 

IN  LIMITS  = 

2 4 

. n n o 

MINIMUM 

VALL'E  OF 

ALL  DATA 

= 

-14 

.800 

MAXIMUM 

VALUE  OF 

ALL  P A T A 

= 

13 

.200 

MOLE  OF 

SATA  = 

- 3 . 

?5  n 

ARITHMETIC  MEAN  r 

- ? . 

553 

STL'.  0 F V I A T TO  M = 

7 . 

334 

VARIANCE 

s 

6 1 . 

376 

V A K I A N C £ 

/ MCA  N = 

* ? 3 • 

ooi 

ESTIMATES  S A M P L F 

MOMENTS 

A ROUT 

THE 

MEAN 

3 

1 6 H . 3 

4 

85  0 3 . 

A3  0.374 


I F 

A 3 

I s 

equal  to 

7ERO,  THF  PIS 

TRIHUTION  IS  SYMMETRICAL 

A 4 

- n 

.540 

I r 

■A  4 

I S 

NEGATIVE 

I)  I STR  I HUT  I ON 

HAS  A FLAT  TOP 

I F 

A 4 

I s 

ZERO,  n I S T R I P U T I n H is 

ALMOST  NOPMAL 

I F 

A 4 

I s 

POSITIVE. 

DISTRIBUTION 

HAS  A HIGH  PEAK. 

I f 

A MY 

p 

S 0 0 A ‘I  I L I T Y 

FUNCTIONS  PO 

NOT  FIT  DATA 

SHOO  FST  THAT  vptj  TRY  THE  FOLLOWING  VALUES 
A M e l = -14.3  0 no 

V|TTM  s 4.1637 

NCELLS=  7 

Figure  C.8.  Statistics  Computations  for  TCAST 
Forecasts  in  Lead  Time  Four 
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This  appendix  tabulates  (see  Tables  D.l,  D.2,  and 
D.3)  tha  following  information  used  in,  and  resulting  from, 
the  statistical  hypotheses  test: 

1.  The  mean  forecast  error  for  each  error  dis- 
tribution. 

2.  The  variance  of  each  distribution. 

3.  The  standard  deviation  for  difference  between 
the  means  for  each  lead  time. 

4.  The  standard  error  of  the  difference  between 
the  means  for  each  lead  time. 

5.  The  number  of  observations  for  each  distri- 
bution. 

6.  The  Student  t Test  Statistic  computed  as  follows: 
a.  For  the  D041  and  TCAST  error  distribution; 


t _ distribution  means  - hypothesized  mean 
distribution  standard  deviation 


b.  For  the  comparison  of  the  two  methods; 


t = 


difference  between  error  distribution  mean 
standard  error  of  the  difference  between 
the  means 


7.  The  confidence  level  at  which  the  null  hypo- 
thesis could  be  rejected. 
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